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The details of Super-Kamiokande-I's solar neutrino analysis are given. Solar neutrino measure- 
ment in Super-Kamiokande is a high statistics collection of *B solar neutrinos via neutrino-electron 
scattering. The analysis method and results of the 1496 day data sample are presented. The final 
oscillation results for the data are also presented. 



I. INTRODUCTION by electron scattering. Due to its unprecedented fiducial 

Super-Kamiokande [Super-K, SK] is an imaging wa- 
ter Cherenkov detector, which detects ^B solar neutrinos 



2 



size of 22.5 kilotons, Super-K has the advantage of mak- 
ing the current highest statistics measurements of solar 
neutrinos. It enables us to determine, with high pre- 
cision, measurements of the solar neutrino flux, energy 
spectrum, and possible time variations of the flux. 

Super-K started taking data in April 1996, and the so- 
lar neutrino results of the first phase of SK, which ended 
in July 2001 and is henceforth referred to as "SK-I," are 
described in this paper. In Section UTI and Section llTTl de- 
tails of the SK detector and its simulation are given. Af- 
ter the event reconstruction method, detector calibration, 
and sources of background are described in Sections IIVI 
Ivl and IVII respectivelv. the data analysis method and re- 
sults are described in Sectior lVlll and IVIIII Finally, in 
Section Hxl the solar neutrino oscillation analysis is dis- 
cussed. 



II. SUPER-KAMIOKANDE DETECTOR 

A. Detector outline 

As has been discussed in much greater detail else- 
where W\, the Super-Kamiokande detector consists of 
about 50000 tons of ultra-pure water in a stainless steel 
cylindrical water tank with 11146 20-inch photomulti- 
plier tubes [PMT's] in the inner detector [ID] and 1885 8- 
inch PMT's in the outer detector [OD]. The diameter and 
height of the SK tank are 39.3 meters and 41.4 meters, 
respectively. The coordinates of the SK tank are defined 
in Figure n In the inner detector, the active photode- 
tector coverage is 40.4% while the remainder is covered 
with black, polyethylene terephthalate sheets, generally 
referred to simply as "black sheet." Signals from PMT's 
are sent through an electronics chain which can measure 
both the arrival time of Cherenkov photons as well as 
the amount of charge they liberate from the phototubes' 
photocathodcs. 
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FIG. 1: Coordinates of the Super-Kamiokande detector. The 
Z-axis is defined as the upward direction, pointing away from 
the center of the earth. 



The Super-K detector is located 1000 meters under- 



ground (2700 meters of water equivalent) in Kamioka 
Observatory, deep within the Kamioka mine in Gifu Pre- 
fecture, Japan. The Observatory is owned and operated 
by the Institute for Cosmic Ray Research [ICRR] , a divi- 
sion of the University of Tokyo. The detector's latitude 
and longitude are 36° 25' N and 137° 18' E, respectively. 
Compared to ground level, the intensity of muons is re- 
duced by about 10~^ at the depth of the SK detector, 
and yielding a downward-going muon rate through the 
detector of about 2 Hz. 

As radon concentrations in the Kamioka mine air can 
exceed 3000 Bq/m"^ during the summer season, there are 
air-tight doors between the Super-K detector area and 
the mine tunnel. The excavated, domed area above the 
cylindrical water tank, called the "SK dome," is coated 
with radon-resistant plastic sheets to prevent radon in 
the surrounding rock from entering the air above the de- 
tector. Fresh air from outside the mine is continuously 
pumped into the SK dome area at the rate of 5 ~ 12 
m^/minute. As a result, the typical radon concentration 
in the SK dome air is 20'--^30 mBq/m'^. 

As we will see in Section IVTl radon can lead to back- 
ground events in the solar neutrino data set. In order 
to keep radon out of the detector itself, the SK tank 
is tightly sealed. Radon-reduced air, produced by a 
special air purification system in the mine, is contin- 
uously pumped into the space above the water surface 
inside the SK tank, maintaining positive pressure. The 
radon concentration of this radon-reduced air is less than 
3 mBq/m'^. 

Finally, the purified water in the SK tank is continu- 
ously circulated through the water purification system in 
the mine at the rate of about 35 tons/hour. This means 
that the entire 50 kiloton water volume of the detector 
is passed through the filtration system once every two 
months or so. 



B. Photomultiplier tubes 

The PMT's used in the inner part of SK are the 20- 
inch diameter PMT's developed by Hamamatsu Pho- 
tonics K.K. in cooperation with members of the origi- 
nal Kamiokande experiment ■ Detailed descriptions of 
these PMT's, including their quantum efficiencies [Q.E.], 
single photo-electron distributions, timing resolutions, 
and so on, may be found elsewhere 0. 

Good timing resolution for light arrival at each PMT 
is essential for event reconstruction. It is around 3 usee 
for single photo-electron light levels. 

The dark noise rate of the PMT's is measured to be 
around 3.5 kHz on average and it was stable over the 
SK~I data taking period as shown in Figure |2Ia). The 
number of accidental hits caused by this dark noise is 
estimated to be about 2 hits in any 50 nsec time window. 
For the solar neutrino analysis, this is corrected during 
the energy calculation as described later in this paper. 

Over time some of the PMT's malfunctioned by pro- 
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ducing anomalously high dark noise rates or emitting 
spark-generated light (a tube which makes its own light 
is called a "flasher"). The high voltage supplied to 
these malfunctioning PMT's was turned off shortly af- 
ter the malfunctions arose, rendering the tubes in ques- 
tion "dead." The number of dead PMT's is shown in 
Figure |2Ib). 

The dark noise rate, plus the numbers of excessively 
noisy and dead PMT's are taken into account in our 
Monte Carlo [MC] detector simulation and also corrected 
for during event reconstruction. 
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FIG. 2: (a) The average dark noise rate of the PMT's used 
in SK-I. The dashed lines show the acceptable range for use 
in the solar neutrino analysis, (b) Number of dead PMT's in 
SK-I. Note that sometimes repairs were possible, usually in- 
volving the replacement of broken high voltage supplies, lead- 
ing to sudden drops in the number of dead PMT's. 



energy [HE] trigger, were set in the ID. The threshold 
number of PMT's for LE and HE triggers were about 29 
and 33 PMTs, respectively. This LE hardware trigger 
threshold corresponded to a 50% triggering efficiency at 
5.7 MeV. In order to avoid hardware trigger efficiency 
issues, the LE analysis threshold was ultimately set to 
6.5 MeV, where the LE hardware trigger efficiency is es- 
sentially 100%. 

From May of 1997, a third ID hardware trigger thresh- 
old called the super low energy [SLE] trigger was added. 
Originally set at 24 PMT's, which provided 50% trig- 
gering efficiency at 4.6 MeV, the addition of the SLE 
trigger served to increase the raw trigger rate from 10 Hz 
to 120 Hz. These SLE events were then passed to an 
online fast vertex fitter, whose operation is described in 
Section HVH 

Since almost of these very low energy events are caused 
by 7S from the rock surrounding the detector and ra- 
dioactive decay in the PMT glass itself, the fast vertex 
fitter was used to reject SLE events with event vertices 
outside the nominal 22.5 kton fiducial volume. This soft- 
ware filtering procedure and its associated online com- 
puter hardware was called the Intelligent Trigger [IT], 
and it filtered the 110 Hz of SLE triggered events down to 
just 5 Hz of SLE events whose vertices fell within Super- 
K's fiducial volume. Thus a total event rate of 14.6 Hz 
was transmitted out of the Kamioka mine for eventual 
offline reduction and analysis. 

In 1999, and again in 2000, the IT system was up- 
graded with additional CPU's. By the end of SK-I it 
provided 100% triggering efficiency at 4.5 MeV, and 97% 
efficiency at 4.0 MeV. Table |2 shows the history of the 
trigger as a function of time. 

The trigger efficiency was checked with both ^^N cal- 
ibration data from our DT generator and events from a 
Ni(n,7)Ni gamma source. It was also continuously moni- 
tored using prescaled samples of real, unfiltered SLE trig- 
gered events. The left plot of Figure |3| shows the typi- 
cal trigger efficiency curve as a function of reconstructed 
energy, which is obtained from the number of hit PMT's 
after various corrections are applied (see Section 11 V CI for 
more details). A trigger efficiency curve as a function of 
true electron energy is calculated by a Monte Carlo sim- 
ulation - it is shown in the right plot of Figure 13 This 
plot is for events whose vertices fall within the fiducial 
volume of the detector (i.e., their true vertex positions 
are more than 2 meters from the PMT wall). 



C. Trigger efficiency 

A data acquisition [DAQ] trigger is generated when- 
ever a certain number of PMT's are fired within a sliding 
200 nsec window. Timing and charge information of each 
fired PMT is digitized once such a hardware trigger is is- 
sued. The ID and OD trigger logic are independent of 
each other. 

When SK-I began data taking in 1996, two threshold 
levels, called the low energy [LE] trigger and the high 



III. SIMULATION 

In the simulation of solar neutrino events in Super- 
Kamiokande-I there are several steps: generate solar neu- 
trinos, determine recoil electron kinematics, generate and 
track Cherenkov light in water, and simulate response of 
electronics. 

In order to generate the ®B solar neutrino spectrum 
in Monte Carlo, the calculated spectrum based on the 
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TABLE I: History of the trigger. The units for rate and live- 
time are Hz and %, respectively. 
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FIG. 3: Trigger efficiency as a function of energy. The left 
plot shows the efficiency as a function of reconstructed energy 
(see Section FlV Cl for the details of reconstructed energy) . The 
May 1997 LE triggers (black circles with black line) and SLE 
triggers (white circles with red line) were calculated using a 
Ni(n,7)Ni gamma source, while the SLE triggers on Septem- 
ber 1999 (triangles with green line) and September 2000 (stars 
with blue line) were calculated using events from a DT 
generator. The right plot shows the efficiency as a function of 
true electron total energy obtained by a Monte Carlo simula- 
tion. Identically colored lines represent the same calibration 
data samples in both plots. 



SK is 325.6 events per day. MC events are generated 
assuming a rate of 10 recoil electron events per minute 
for the full operation time of SK-I, yielding a total of 
24,273,070 simulated solar neutrino events. 
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FIG. 4: (a) Input distributions of *B (solid) and hep (dashed) 
solar neutrino energies, (b) The cross section of the interac- 
tion for Ve (solid line) and i/^,T (dashed line) with electrons 
as a function of neutrino energy, (c) The spectrum of recoil 
electrons scattered by *B and hep solar neutrinos. 



decay measurement of Ortiz et al. |^ was used. Fig- 
ure ^a) shows the input solar neutrino energy distribu- 
tion for neutrinos. For the uncertainty in the spec- 
trum we have adopted the estimation by Bahcall et al . 

In the next step, the recoil electron energy from the 
following reaction; 

ly + e ^ ly + e (3.1) 

is calculated. Fig^Jb) shows the differential cross section. 
The (j^e,e) scattering cross section is approximately six 
times larger than (i'p.T-,e), because the scattering of ly^ 
on an electron can take place through both charged and 
neutral current interactions, while in case of Vf^ r only 
neutral current interactions take place. The radiative 
corrections in the scattering are also considered 6] . Fig- 
ure ^c) shows the expected spectrum of recoil electrons 
in SK without neutrino oscillations. 

Assuming the ^B flux of BP2004 (5.79 x 10^ / cm'^ / sec) 
and before taking into account the detector's trigger re- 
quirements, the expected number of j/g scatter events in 



We have used GEANT3.21 for simulation of particle 
tracking through the detector. Since the tracking of 
Cherenkov light is especially important in the Super-K 
detector, parameters related to photon tracking are fine 
tuned through the use of several calibration sources. 

Figure |S1 shows the wavelength dependence of various 
water coefficients in our MC. Rayleigh scattering is dom- 
inant at short wavelengths with a 1/A^ dependence. The 
A dependence of absorption and Mie scattering are em- 
pirically set to 1/A^ at shorter wavelengths, while the 
absorption for longer wavelengths are taken from a sepa- 
rate study 7] . The absorption and scattering coefficients 
in MC are tuned using LINAC calibration data (see Sec- 
tion ^^2J so as to match the MC and data energy scale 
in each position in the detector. 

The water quality in SK changes as a function of time 
as shown in Figure [T7I this change in water quality was 
taken into account in the MC simulation. By compar- 
ing photon arrival timing distributions using calibration 
data, it was found that the change in water attenuation 
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length is mainly due to change in the absorption coeffi- 
cient. So, we fixed Rayleigh scattering and Mie scattering 
coefficients to be constant over the entire data taking pe- 
riod and vary only the absorption coefficient in the MC 
simulation. Figure shows the result. The solid line 
shows the translation of water transparency into the ab- 
sorption coefficient. The lower panel shows the deviation 
of the peak of the reconstructed energy from the peak of 
the input energy as a function of water transparency — 
the deviation is less than ±0.2%. The coefficients of each 
process at shorter wavelength are summarized as follows: 



8.00 X 107A''[nTO] (1/m) 
1.00 X 10VA'*[nTO] (1/m) 
(2.74 - 9.27) X 10VA^[7ito] (1/m) 
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Mie scat., 
absorption. 
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FIG. 5: Wavelength dependence of the water parameter coef- 
ficients: absorption (solid), Rayleigh scattering (dashed) and 
Mie scattering (dotted). The absorption coefficient is also a 
function of water transparency. The filled region shows the 
range of this parameter as water transparency is changed, 
where the two solid lines define the SK-I minimum (73 me- 
ter) and maximum (98 meter) values. 



IV. EVENT RECONSTRUCTION METHOD 
A. Vertex 
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FIG. 6: Translation of the water transparency parameter into 
a short wavelength absorption coefficient. Absorption is taken 
to be proportional to A"* for these wavelengths. Light scat- 
tering is assumed to be constant. The lower panel shows the 
energy deviation as a function of input momentum. The inner 
error bar of each point is the statistical uncertainty, while the 
outer error bar represents the spread of the six samples using 
5, 6, 8, 10, 12, and 15 MeV/c as the input momenta. 
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FIG. 7: Search grid for the vertex reconstruction in the xy 
plane and in z. The shaded area is the SK inner detector 
volume. The dashed line indicates the fiducial volume. The 
white stars (white lines) are the grid points used by the stan- 
dard vertex reconstruction. The black circles are used by the 
second vertex reconstruction. 



Electrons in the energy region of interest for solar neu- 
trinos (below 20 MeV) can travel only a few centime- 
ters in water, so their Cherenkov light is approximately 
a point source. The reconstruction of this vertex relies 
solely on the relative timing of the "hits," i.e., PMT's 
struck by one or more Cherenkov photons. Since the 
number of observed Cherenkov photons and therefore the 
likelihood of a multiply hit PMT is comparatively small, 
about seven recorded photons per MeV of deposited en- 
ergy, the pulse heights of the hits typically follows a 
one photo-electron distribution and yields no information 
about the light intensity nor the distance to the source. 



For the same reason it is also impossible to separate re- 
flected or scattered Cherenkov photons and PMT dark 
noise from direct light based upon PMT pulse height. 
The vertex reconstruction assumes that a photon origi- 
nating from vertex v and ending in hit hi is traveling on 
a straight path and therefore takes the time \v — hi\/c 
where c is the group velocity of light in water (about 
21.6 cm/ns). Therefore, assuming only direct light, the 

effective hit times — ti — \v — hi\/c at vertex v of all 
PMT's should peak around the time of the event with 
the width of the PMT timing resolution for single photo- 
electrons (about 3 ns). Light scattering and reflection (as 
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well as dark noise, pre-, and afterpulsing of the PMT's) 
introduces tails in the distribution which will strongly 
bias the reconstruction if a is used to evaluate the 
goodness of fit g for a given vertex. Therefore, we use a 
"truncated x^" 

(4.1) 

1=1 

where to is the center of a 10 ns wide search window in r 
which maximizes the number of hits inside. The vertex 
is reconstructed by the maximization of this goodness 
varying v. 

We use three different vertex reconstruction algo- 
rithms. Our standard vertex fit is used to reconstruct 
event direction and energy and to compute the likelihood 
of an event to be due to spallation. Only events inside 
the fiducial volume (2 m away from the closest PMT) are 
considered. The Super-K background rate below about 
7 MeV rises rapidly with falling energy; most of this back- 
ground is due to light emanating from or near the PMT's 
themselves and is reconstructed outside the fiducial vol- 
ume. However, the rate of misreconstructed events (back- 
ground events outside the fiducial volume which get re- 
constructed inside) increases rapidly with falling energy 
due to long resolution tails. We therefore reconstruct 
the vertex using a second vertex fitter whose vertex dis- 
tribution has smaller tails and then accept only events 
reconstructed inside the fiducial volume by both fits. To 
reduce the rate of SK events stored to tape, we remove 
SLE events which reconstruct outside of the fiducial vol- 
ume online. Unfortunately, both offline fitters are too 
slow to keep up with the trigger rate (see Section Hi C|l 
of SLE events and so we are forced to use a fast recon- 
struction to pre-filter these events. If the vertex of the 
fast online fit is inside the fiducial volume the event is 
reconstructed later by the other two fits. 

To ensure convergence of the maximum search, the 
standard fit first evaluates the vertex goodness on a 
Cartesian grid (see Figure [7)). For a reasonably speedy 
search in spite of Super-Kamiokande's large size, a coarse 
grid of 57 points with a grid constant of 397.5cm is cho- 
sen for the xy plane with nine such layers in z separated 
by 380 cm. Since the grid is coarse, the timing resolution 
at is artificially set to 10 ns to smear out the maximum 
goodness. After the single coarse grid point with the 
largest goodness is identified, the goodness is calculated 
at the 27 points of a 3 x 3 x 3 cube centered on this grid 
point with a timing resolution of <7t = 5 ns and a spa- 
tial grid separation of 156 cm. If the largest resulting 
goodness is not in the center, the cube is shifted so that 
the largest goodness is in the center of the new cube. 
Otherwise, the cube is contracted by a factor of 2.7. If 
the spatial grid separation of the cube reaches 5 cm the 
maximum search is finished and the reconstructed vertex 
is the point with the largest found goodness. 

To limit the bias of the vertex reconstruction due to 
the tails in the timing residuals ti — t^, the hits that go 
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FIG. 8: Timing constraint for pairs of hits. From the triangle 
relation cti + /S.x > ctj = cti + Actij follows At^j < Aa;/c = 
\h, - hj\/c. 



in the vertex fit are selected. Since the event vertex is 
not known, the hit selection must be based on the abso- 
lute hit times ti (and the spatial distribution of the hits 
hi). The hit selection of the standard reconstruction first 
finds the start time of a sliding 200 ns window (the time 
required for a photon to cross the diagonal of the detec- 
tor) which contains the largest number of hits. The rate 
of late/early hits (background) per ns is determined by 
counting the hits outside the window. From this rate, 
the background inside the window is estimated. The size 
of the window is then reduced in an attempt to optimize 
the direct light signal divided by the square root of the 
late/early hit background. The hits in the resulting, op- 
timal timing window are those selected to compute the 
vertex goodness. 

The second vertex reconstruction (to reduce misrecon- 
structed events originating from the PMT's) modifies the 
hit selection and the search grid. For two direct light 
hits i and j, the timing difference Ai^ is limited to 
A^ij < \hi — hj\/c (see Figure ISJ. We select the largest 
set of hits whose hit pairs obey Aty < \hi — hj\/ c after 
eliminating "isolated hits" (hits which are further away 
than 1250 cm or further away than 35 ns from the near- 
est neighboring hit). With these selected hits we perform 
a grid search with a circular grid of 60 points on the xy 
plane (grid constant 397.5 cm) and nine such planes in z 
separated by 380 cm with a timing resolution of 9.35 ns. 
The best fit point is interpolated from the grid points 
with the largest goodness. After that, 27 points of a 
cylindrical section around that point with an initial grid 
constant of 147.1 cm and 5 ns resolution are tested. As 
in the standard reconstruction, the section is moved (if 
the center point doesn't have the largest goodness) or 
reduced in size by a factor of 0.37 (if it does), and the 
search is finished once the grid constant falls below 5 cm. 

The fast fit, used to pre-filter SLE events online, also 
eliminates isolated hits (i.e., all other hits separated by 
either more than 10 m or more than 33.3 ns) to reduce 
the effects of dark noise and reflected or scattered light. 
Then the absolute peak time is estimated by maximiz- 
ing the number of PMT hits within a 16.7 ns wide sliding 
timing window. PMTs that are in the interval (—33.3 ns. 
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FIG. 9: Search grid for the fast fit. 18 points around the 
input vertex are tested, the orientation is fixed. The grid size 
is adjusted between one and eight meters (see Table . 



Step Size [cm] 


to ± at [ns] 


do ± 


twin [ns] 


1 


800 


2.1 ± 10.2 


0.534 ±0.134 


(-66.7,100.0) 


1 




2.7 ±6.1 


0.534 ±0.134 


(-40.0, 66.7) 


2 


400 


2.3 ±6.0 


0.534 ± 0.134 


(-33.3, 50.0) 


3-5 




2.3 ±6.0 


0.534 ±0.134 


(-25.0, 33.3) 


1 




1.5 ±4.7 


O.lg for a < 0.2, a > 0.9 


(-20.0,33.3) 


2 


200 


2.2 ±4.5 


g for 0.3 < a < 0.85 


(-20.0,33.3) 


3-5 




1.1 ±4.1 


0.5g otherwise 


(-20.0,26.7) 


1-8 


100 


0.6 ±3.2 


none 


(-13.3,16.7) 



TABLE 11: Search parameters for the fast fit. 



100 ns) with respect to the peak time are selected. An 
initial vertex is calculated by shifting a simple average 
of the selected PMT positions 2 m toward the detector's 
center. The time at the event vertex is also determined; 
this time must be smaller than the absolute peak time 
yet not differ from it by more than 117 ns. The initial 
vertex time is chosen to be 58.3 ns before the absolute 
peak time. It is then corrected to the average of all time- 
of-flight subtracted PMT times that are in the interval 
(—133.3 ns, 200 ns) around the initial time. The magni- 
tude of the anisotropy 

^selected liits 

is also calculated at this stage. Next, 18 points around 
the initial vertex are tested with an initial search radius 
of 8 m (see Table EJ, while the goodness g is modified to 
take into account the anisotropy. At first, the goodness 
is multiplied by the factor e~''-^'('^~'^o)/o-a) _ Later in the 
search, this factor is replaced by a table. The times are 
shifted by the expected mismatch between the average 
vertex time and the vertex peak time Iq. Only hits in- 
side a time interval around the shifted vertex time are 
considered. The search radii, time shifts, time resolu- 
tions, anisotropy factors, and time intervals are listed in 
Table HH 

As shown in Figures ITUl and 1771 the performance of all 
three fitters is evaluated using Monte Carlo. The dis- 
tribution of the distance between the reconstructed and 
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FIG. 10: Distance from reconstructed to correct vertex for 
*B Monte Carlo events with total recoil electron energy above 
4.5 MeV. The inserted panel magnifies the distributions near 
zero on a hnear scale. 
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FIG. 11: Vertex resolution (68.2% of reconstructed events 
reconstruct inside a sphere of radius a from correct vertex) of 
*B Monte Carlo events as a function of total recoil electron 
energy. The curves on the right-hand side show the fraction 
of reconstructed events as a function of total recoil electron 
energy. 



correct vertex is analyzed for each fit (see Figure^)); 
the vertex resolution is extracted as the distance which 
contains 68.2% of all reconstructed vertices. For recoil 
electrons above 4.5 MeV, the fast fit shows the worst res- 
olution (115 cm). However, it has the lowest rate of very 
distant misfits (> 30 m). The resolution of the standard 
fit is 102 cm, while the second fit's resolution is 94 cm. 
Figure ^2 shows the vertex resolution as a function of 
generated total recoil electron energy. Since the vertex 
reconstruction occasionally fails (due to a number of se- 
lected hits that is too small) , the efficiency to reconstruct 
the event anywhere at all is also plotted. 
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B. Direction 

Since the recoil electron preserves the direction of the 
solar neutrino, directional reconstruction is important for 
solar neutrino analysis. This characteristic directionality 
is used to extract the solar neutrino signal in Super-K. 
To calculate the direction, a maximum likelihood method 
using the Cherenkov ring pattern is adopted. The likeli- 
hood function is 



L{d) EE ^log(/(cos6'd^r))j X 



cos 9i 



(4.3) 



where iVso is the number of hit PMT's with residual times 
within a 30 nsec window and f{cos9dir) is the function 
that represents the distribution of the opening angle be- 
tween the particle direction and the vector from the re- 
constructed vertex to the hit PMT position made by MC. 
A plot of f{cos6dir) for 10 MeV electrons is shown in 
Figure El The distribution is broad with the peak at 
42° because of the effects of electron multiple scatter- 
ing and Cherenkov light scattering in water. 9i is the 
opening angle between the direction of the vector from 
the reconstructed vertex to the i-th hit PMT position 
and the direction that PMT is facing. a{9i) is the ac- 
ceptance of PMT's as a function of 6i which is made by 
MC. The direction is reconstructed as Eq. 14.3|l becomes 
maximum using a grid search method whose step sizes 
are 20°, 9°, 4°, 1.6°. 




cose.: 



FIG. 12: The distribution of the opening angle between the 
direction of the generated particle and the vector from the 
reconstructed vertex to the hit PMT position. This plot is 
made for 10 MeV electron MC events. 

The quality of the directional reconstruction is esti- 
mated from the differences between the generated and 
reconstructed directions using MC. Figure [TSI shows the 
directional resolution dependence on energy within the 
fiducial volume; it uses directions which are generated 
uniformly. Here, the resolution is defined as the 68% 
point for the angle difference distribution. The angular 
resolution for 10 MeV electrons is about 25 degrees. 
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FIG. 13: The directional resolution's dependence on energy, 
determined using Monte Carlo events. 



C. Energy 

The energy of a fully-contained charged particle in 
Super-Kamiokande is approximately proportional to the 
number of generated Cherenkov photons, and thus is also 
proportional to the total number of photo-electrons in the 
resulting hit PMT's. When the particle's energy is low it 
is also proportional to number of hit PMT's because in 
such a case the number of Cherenkov photons collected 
by any given PMT is almost always zero or one. In or- 
der to avoid the effect of noise hits with higher charge, 
the number of hit PMT's {Nhu) with some corrections 
is used for energy determination for the solar neutrino 
analysis. In order to reject accidental hits due to dark 
noise in the PMT's, only the hit PMT's with residual 
times within a 50 nsec window are used for calculating 
Nhit- Moreover, we applied several corrections to Nhu, 
yielding an effective number of hits {N^f / ) which has the 
same value at every position in the detector for a given 
particle energy. These corrections account for the varia- 
tion of the water transparency, the geometric acceptance 
of each hit PMT, the number of bad PMT's, the PMT 
dark noise rate, and so on. The equation for N^ff is: 



^^ff ~ ^ {{-^i + ^tail — f-dark) X 



i=l 

-Rcover 



N, 



normal 



X exp 



A (run) 



X G,(t)} (4.4) 



where Xi is an occupancy used to estimate the effect of 
multiple photo-electrons, etaii is the correction for late 
hits outside the 50 ns window, and Cdark is for dark noise 
correction. The definition of X,- is as follows: 



Xi 



Xi <i 
Xi = I 



(4.5) 
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where Xi is the ratio of the number of hit PMT's to the 
total number of PMT's in a 3 x 3 patch around the i-th hit 
PMT. This correction estimates the number of photons 
which arrived at the i-th hit PMT by using the number 
of hit PMT's surrounding it. Here, the ratio of the unhit 
PMT's in this patch of nine tubes is 1 ~ Xi. 

The second factor is the bad PMT correction, where 
Nail is the total number of PMT's, 11146, and Nnormai is 
the number of properly operating PMT's for the relevant 
subrun. 

The third factor is the effective photo coverage. The 
average of the photo coverage, which is the ratio of the 
area covered by PMT's to all inner detector area, is 
Rcavei — 0.4041 . However, the effective photo cover- 
age changes with the incident angle of the photon to the 
PMT. Therefore, we applied a coverage correction, where 
S{9i,(j)i) is the photo coverage from the directional van- 
tage point of Oi,(l)i. Figure ITU shows this function. 

The fourth factor is the water transparency correction, 
where is the distance from the reconstructed vertex to 
the i-th hit PMT position, and A is water transparency. 

The fifth factor, G{t)i, is the gain correction at the 
single photo-electron level as a function of the time of 
the manufacture of each PMT. 




FIG. 14: Function of the effective photo-coverage's depen- 
dence on the incident angle to a PMT. 



Finally, the conversion function from N^f f to energy is 
determined by uniformly generated Monte Carlo vertices. 
This is total energy; all mention of "energy" hereafter 
refers to total energy, i.e. including the scattered elec- 
tron's rest mass and momentum. The relation to energy 
is calculated by the following equation; 

E = a + /5iVe// ( 1 - 7A^e// ( 1 - ^ A^e// ( 1 - eA^e// ) ) ) , (4.6) 

here a = 0.735, /? = 0.134, 7 = 6.049 x 10"^, 5 = 6.441 x 
lO^'^, and e — 1.541 x 10~^. The typical conversion factor 
from Neff to energy is 6.97. The energy resolution is also 
estimated by the same MC described above. Figure [T^ 
shows the energy resolution as a function of energy; it is 
14.2% for 10 MeV electrons. 




5 10 15 20 25 

Energy (MeV) 

FIG. 15: The energy resolution's dependence on energy as 
determined using Monte Carlo events. 



D. Muon 

Precise reconstruction of cosmic ray muon tracks which 
penetrate the Super-K detector is needed for the so- 
lar neutrino analysis. This is because nuclear spallation 
events induced by these cosmic ray muons are a dominant 
background to the solar neutrino signal, and the corre- 
lation in time and space between spallation events and 
their parent muons is extremely useful in rejecting the 
spallation background. Hence, the track position preci- 
sion requirement for our muon fitter is about 70 cm since 
the vertex resolution for low energy events with energies 
around 10 MeV is also about 70 cm. 

The muon fitter has three components; an "initial fit- 
ter," a "TDC fitter," and a "geometric check." 

The initial fitter assumes the position of the first fired 
PMT is the entrance point of a cosmic ray muon, and 
the center of gravity of the saturated tubes' positions 
as the exit point. The PMT at the exit point has the 
most photo-electrons [p.e.'s] of all PMT's; the expected 
value is up to 500 p.e.. However, our front end electronics 
saturate at ~ 230 p.e. in a single tube, and so typically 
several dozen PMT's near the exit point are saturated. 
Therefore, the reconstructed exit point is not the position 
of the most-hit PMT, but rather the center of gravity of 
the saturated tubes' positions. 

For the entrance position search, we consider the dark 
noise of the PMT's. For SK-I the dark noise rate of a 
typical PMT was ^ 3.5 kHz, though it could rise to > 10 
kHz in the case of a "noisy PMT." In the initial fitter, 
the following methods are used to reject dark noise hits: 

1. Charge information of PMT 

The fitter requires that the PMT at the entrance 
point have more than 2 p.e., because typical dark 
hit PMT's have less than 1 p.e. while the amount 
of charge which they receive from a muon is several 
p.e.'s. 
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2. Timing information of eight neighbor PMT's 
Cosmic ray muons emit lots of Cherenkov light 
(~ 340 photons/cm), so several PMT's near the 
entrance point should record photons at the same 
time. The fitter therefore requires that the PMT at 
the entrance point have more than five fired neigh- 
bor PMT's within 5 ns. 

The results of the "initial fitter" are then used as initial 
values by the next stage. 

The TDC fitter is a fitter based on a grid search 
method using timing information. In this fitter, a good- 
ness of fit is also defined to obtain the best track of the 
muon: 

g(v)^2_^e ^\7^) ^ (4.7) 

i=l 

where the definitions are the same as for the "truncated 
X^" in Section FlV Al except for the factor (/), which is 
found to be 1.5 by Monte Carlo. The TDC fitter surveys 
a circle of radius 5.5 m from the exit point obtained by 
initial fitter and then obtains the direction which maxi- 
mizes the goodness. 

The muon fitter determines the muon track by using 
two independent parameters, timing and charge, and the 
two results are sometimes different. For example, in the 
case of muon bundle events, PMT's near the exit point 
fire at early times, because the velocity of a cosmic ray 
muon is faster than that of light in water. As a result 
of these early hits, the exit point can sometimes be mis- 
taken for an entrance point. In order to reject this kind 
of misfit, the geometric check fitter estimates the consis- 
tency of both previous fitters by requiring the following 
after the TDC fit: 

• There is no saturated PMT within 3 m from the 
entrance point. 

• There should be some saturated PMT within 3 m 
from the exit point. 

The muon events which do not satisfy these requirements 
are regarded as failed reconstructions. 

We estimated the performance of the muon fitter by us- 
ing Monte Carlo events. Figure [TBI shows the vertex and 
angular resolutions. For the entrance point, the 1-a ver- 
tex resolution is estimated to be 68 cm. This value corre- 
sponds closely to the spacing interval between each PMT, 
70.7 cm. The vertex resolution of exit point is 40 cm. 
Therefore, the distance between the actual muon track 
and the fitted one is estimated to be less than 70 cm, 
which means that this muon fitter satisfies our precision 
requirement. The angular resolution is estimated to be 
1.6°. 
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FIG. 16: Performance of the muon fitter. Left (right) figure 
shows the vertex (angular) correlation between generated and 
reconstructed Monte Carlo muon tracks. The area to the 
left of the dashed line in both plots shows the 1-cr (= 68%) 
included region. 

V. CALIBRATION 
A. Water transparancy measurement 

Cherenkov photons can travel up to 60 m before reach- 
ing PMT's in Super-K, and so light attenuation and scat- 
tering in water directly affects the number of photons 
that are detected by the PMT's. Since the energy of 
an event is mainly determined from the number of hit 
PMT's, the water transparency [WT] must be precisely 
determined for an accurate energy measurement. 

The water transparency in SK is monitored continu- 
ously by using the decay electrons (and positrons) from 
cosmic ray /i events that stop in the detector: /i^ — s- 
e~ + + v^, and ^ e+ + i^e + S'^t- At SK's under- 
ground depth (2700 m.w.e.), cosmic ray muons reach the 
detector at a rate of ~ 2 Hz. Approximately 6000 /i's 
per day stop in the inner detector and produce a decay 
electron (or positron) . In order to monitor the WT effec- 
tively, it is important to have a pure sample of /i-e decay 
events. Several criteria are applied to select these events: 

• The time difference [AT] between the stopping /i 
event and the /^-e decay candidate must satisfy: 
2.0 /isec < AT < 8.0 ^sec. 

• The reconstructed vertex of the /i-e decay candi- 
date must be contained within the 22.5 kton fidu- 
cial volume. 

• N}iit must be at least 50. 

These criteria select ~1500 /i-e decay events daily, which 
is sufficient statistics to search for variation in the WT 
in one week. The average energy of the /i-e decay events 
is ~ 37 MeV, which is much higher than that of solar 
neutrinos (< 20 MeV). Therefore, the /i-e decays cannot 
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be used for the absolute energy calibration. However, 
they can be used to monitor the stability of water trans- 
parency and energy scale over time. 

In order to remove the effects of scattered and reflected 
light, hit PMT's are selected by the following criteria: 
(1) PMT's must have timing that falls within the 50 nsec 
timing window, after the time-of-flight [TOF] subtrac- 
tion, (2) PMT's must be within a cone of opening angle 
32° ~ 52° with respect to the direction of the fi-e decay 
event. A plot of the number of hit PMT's vs. distance 
from the vertex to the PMT is made using the selected 
PMT's, fit with a linear function, and the inverse of the 
slope gives the water transparency. 

Figure [TTT a^ shows variations in the WT as a function 
of time. Each point on the plot represents WT for one 
week. In order to reduce the effects of statistical fluctu- 
ations on the weekly measurement, the WT for a given 
week is defined as a running average over five weeks of 
data. Figure [rTr bl plots the mean N^fj value for ^-e de- 
cay events in SK-I. From this figure, it can be seen that 
the energy scale has remained stable to within ± 0.5% 
during the SK-I runtime. 
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FIG. 17: (a) Time variation of the measured water trans- 
parency during SK-I. Stability of the SK-I energy scale as a 
function of time, (b) Each point represents the mean Neff 
value for 28 days of n-e data. The solid line represents the av- 
erage Neff for all ^-e events in SK-I; the dashed and dotted 
lined represent, respectively, ± 0.5% and ± 1.0% deviations 
from the average. 



B. Energy calibration 

The energy scale calibration is necessary to correctly 
convert the effective number of hits {Neff) to the total 
energy of the recoil electrons induced by solar neutrinos. 
LINAC and DT calibration are discussed in this section. 



1. LINAC 

The primary instrument for energy calibration is an 
electron linear accelerator [LINAC] . The LINAC calibra- 
tion of Super-K has been discussed in detail elsewhere |3] ■ 
The LINAC is used to inject downward-going electrons 
of known energy and position into the SK tank. The 
momentum of these electrons is tunable, with a range 
of 5.08 MeV to 16.31 MeV; this corresponds well with 
the energies of interest to solar neutrino studies. LINAC 
data is collected at nine different positions in the SK 
tank, as shown in Figure IT51 The data are compared to 
MC, and MC is adjusted until its absolute energy scale 
agrees well with data in all the positions and momenta. 
Once adjusted, this MC is extrapolated to cover events 
in all directions throughout the entire detector. 

The electrons are introduced into SK via a beam pipe 
whose endcap's exit window is a 100 /im-thick sheet of 
titanium. This allows electrons to pass through with- 
out significant momentum loss but prevents water from 
entering the beam pipe. 
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FIG. 18: The LINAC system at SK. The dotted line repre- 
sents the fiducial volume of the detector and the black dots 
indicate the positions where LINAC data were taken. 



Seven different momentum data sets are taken at each 
LINAC position. The absolute energy in the LINAC sys- 
tem is measured with a germanium [Ge] detector. Table 
^3 shows the different momenta used in the LINAC cali- 
bration. The average occupancy at the beam pipe endcap 
is set to about 0.1 electrons/spill. The reason for this low 
occupancy is to reduce the number of spills which include 
multiple electrons. 

The energy scale and resolution obtained by LINAC is 
compared to MC at each position and momentum. Fig- 
ure liyf a') shows the deviation in the energy scale between 
data and MC. FigureE^b) shows the average over all the 
positions from (a). The deviation at each position is less 
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beam momentum 
(MeV/c) 


Ge energy 
(MeV) 


in-tank energy 
(MeV) 


5.08 


4.25 


4.89 


6.03 


5.21 


5.84 


7.00 


6.17 


6.79 


8.86 


8.03 


8.67 


10.99 


10.14 


10.78 


13.65 


12.80 


13.44 


16.31 


15.44 


16.09 



TABLE III: LINAC beam momentum. The second column 
gives the energy measured in the Ge cahbration system. The 
last column lists the total energy of the electrons after leaving 
the beam pipe. 



O 

I 0.05 



1,0.05 


-0.05 







I 

i 


----ii'---- 


-If: 


1 ....... 1 ....... 1 . 

(b) : 


■ * * * 




i 


t 


i 





4 6 8 10 12 14 16 

Energy (MeV) 



o 
< 

z 0.01 


(a) f 




i 1 : 


< 

Z-0.01 




o 

1,0.01 


i i i 1 i i i 1 i i i 1 i i i 1 i i i 1 i i i 1 i 
(b) : 





- i i i 






-0.01 





4 6 8 10 12 14 16 

Energy (MeV) 



FIG. 19: The relative difference between the reconstructed 
energy of LINAC data and the corresponding MC. (a) The 
results are shown for all positions and beam momenta, (b) 
Averaged over all positions. The dashed lines show ± 0.5%. 
See Fig. 1181 for the positions. 



than 1% and the position averaged momentum depen- 
dence of the deviation is less than 0.5%. 

Figure [211 (a) shows the deviation of the energy reso- 
lution at each energy and position in LINAC cahbration 
between data and MC, and (b) shows these deviations av- 
eraged over all positions. The difference in reconstructed 
energy resolution between data and MC is less than 2.5%. 

The absolute energy scale of the detector was tuned by 
using the LINAC calibration data. The uncertainty of 
the absolute energy scale comes from the uncertainties of 
water transparency measurement when the LINAC cali- 
bration data was taken (0.22%), position dependence of 
the energy scale (0.21%), time variation of the energy 
scale obtained by LINAC calibration (0.11%), tuning ac- 
curacy of MC simulation (0.1%), electron beam energy 
determination in LINAC calibration (0.21%), and direc- 
tional dependence of energy scale since the LINAC beam 
is only downward-going (0.5%). Adding those contribu- 
tions in quadrature, the total uncertainty of Super-K-I's 
absolute energy scale is estimated to be 0.64%. 



FIG. 20: The reconstructed energy resolution from LINAC 
data and MC. (a) shows the deviation of the resolution in 
each energy and position in LINAC calibration between data 
and MC. (b) Averaged over all positions. The dashed lines 
show ± 2.5%. See Fig. ll8l for the positions. 



2. ^''N from the DT Generator 

Even though the energy scale is determined from 
LINAC calibration and is fixed, ^^N was also used as 
a calibration source. Since events from ^^N decay are 
isotropic, they are useful in probing the directional de- 
pendence of the energy scale. Also, because of the porta- 
bility of the DT generator we were able to probe the en- 
ergy scale at many positions in the detector on a monthly 
basis. 

With a half-life of 7.13 seconds, the Q-value of the de- 
cay of ^^N is 10.4 MeV, and the most probable decay 
mode produces a 6.1 MeV 7 ray together with a /3 decay 
electron of maximum energy 4.3 MeV. Man-made ^^N 
was obtained using a deuterium-tritium neutron genera- 
tor ["DT generator"]. The DT generator uses the fusion 
reaction '^H +^ H He + n to produced 14.2 MeV 
neutrons. A fraction of the neutrons collide with ^^O to 
produce ^^N. The DT generator thus provided us with a 
large sample of essentially background free ^^N data for 
use in energy scale calibration. More details about the 
DT generator are given elsewhere . 

The ^^N decay energy spectrum is a superposition of 
several 7 ray lines and (3 continua of various end points. 
The reconstructed energy spectrum has a peak around 
Epeak — 6.9 ^ 7.0 MeV. The shape of the peak region 
(5.5 ^ 9.0 MeV) is approximated by a Gaussian with 
a width of 1.6 ~ 1.7 MeV. The deviation of this energy 
peak between data and MC is measured. FigureETlshows 
31 data sets of DT results during two years of running, 
and it shows the deviation is within ±1% when averaged 
over all positions. 

The direction dependence of the energy scale is mea- 
sured as follows. First, the deviation between data and 
MC is calculated for all directions. Then it is calculated 
for two divided data samples: events in the upward or 
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FIG. 21: The deviation between data and MC for the volume- 
averaged DT energy scale as a function of time. The DT 
generator was available for the last two years of SK-I, starting 
in July, 1999. The dashed line shows the ±0.64% of energy 
scale uncertainty, which is estimated by LINAC calibration 
as discussed in Sectior fV B II 
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C. Vertex calibration 

Figure [221 shows the vertex resolution in each energy 
and position in LINAC and Monte Carlo. The differences 
of vertex resolution between data and MC is less than ±5 
cm. 

The vertex shift in the reconstruction is estimated us- 
ing a Ni(n,7)Ni gamma source, Q because the gamma 
ray is emitted in almost uniform directions. The vertex 
shift is defined as a vector from an averaged position of 
the reconstructed vertex of the data to that of a corre- 
sponding MC. Table Hvl shows the vertex shift at several 
source positions. The systematic error for the solar neu- 
trino flux as a result of vertex shift (which could move 
events in or out of the fiducial volume) is evaluated from 
these values, and it is ±1.3%. 
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FIG. 23: The reconstructed vertex resolution from LINAC 
data and MC. (a) shows the deviation of the resolution in 
each energy and position in LINAC calibration between data 
and MC. (b) averaged over all positions. The dashed lines 
show ± 5 cm. See Fig. 1181 for the positions. 
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FIG. 22: Directional dependence of energy scale by DT cal- 
ibration. Upper figure shows Rdown — R as a. function of 
time, lower figure shows the same for Rup — R- Note that 
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D. Direction calibration 



downward directions. It should be noted that in all cases, 
the energy scale is averaged over the detector volume. 
Figure shows plots comparing the energy scale ob- 
tained from using just downward- or upward-going events 
to that determined by using the entire sample. The dif- 
ference is within ±0.6%. 



The angular resolution is calibrated by LINAC data 
and MC. Figure |^ shows the angular resolution vari- 
ation in each energy and position between LINAC and 
Monte Carlo events. The difference in angular resolution 
between LINAC data and MC is less than ±0.5 degree. 

We have applied this difference as a correction factor to 
the expected signal shape, then taken the same amount 
of the correction as our systematic error due to angular 
resolution. This systematic error is 1.2%. 
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Position 


Aa; 


Ay 


A^ 


(35.3,-70.7,-1200) 


-1.8 


-1.9 


-2.8 


(35.3,-70.7,0) 


0.6 


-0.5 


-2.8 


(35.3,-70.7,-1-1200) 


-1.1 


-1.7 


0.6 


(35.3,-70.7,+1600) 


0.0 


-3.2 


-3.5 


(35.3,-1555.4,-1200) 


-5.5 


8.6 


-9.1 


(35.3,-1555.4,0) 


-5.5 


23.4 


-3.0 


(35.3,-1555.4,4-1200) 


-0.5 


7.0 


6.0 


(35.3,-1555.4,4-1600) 


-3.3 


7.6 


2.7 


(35.3,-1201.9,-1200) 


-4.2 


5.4 


-9.1 


(35.3,-1201.9,0) 


-1.4 


15.5 


-2.2 


(35.3,-1201.9,-^1200) 


2.5 


8.9 


4.2 


(1520.0,-70.9,-1200) 


-5.7 


-1.4 


-10.3 


(1520.0,-70.9,0) 


-18.6 


-3.0 


-4.0 


(1520.0,-70.9,+1200) 


-12.6 


-1.0 


6.2 


(-35.3,-1555.4,-1200) 


3.7 


-11.4 


-8.6 


(-35.3,-1555.4,0) 


6.6 


-19.5 


-2.0 


(-35.3,-1555.4,+1200) 


1.6 


-20.3 


7.2 



TABLE IV: The vertex shift measured by Ni(n,7)Ni gamma 
source calibration. The unit is cm. 
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FIG. 24: The reconstructed angular resolution from LINAC 
data and MC. (a) shows the deviation of the resolution in 
each energy and position in LINAC calibration between data 
and MC. (b) is averaged over all positions. The dashed lines 
show ± 0.5 degree. See Fig. 1181 for the positions. 



VI. BACKGROUND 
A. Low energy backgrounds 

The main background sources below about 6.5 MeV 
for the solar neutrino events are (1) events coming from 
outside fiducial volume, and (2) ^^^Rn. 

Figure 1231 shows a typical vertex distribution of the 
low energy events before the Intelligent Trigger selection. 
Most of the low energy events occur near the inner de- 
tector's wall region. They originate from radioactivity 
of the PMT's, black sheet, PMT support structures, and 
mine rocks surrounding the SK detector. Though the 
fiducial volume is selected to be 2 m from this inner de- 




3000^ 
X 10 ' 



(cm^) 



FIG. 25: Typical vertex distribution of the low energy events 
before the Intelligent Trigger selection. The analysis energy 
threshold for this plot is 5.0 MeV. The dashed line shows the 
fiducial volume edge. 



tector wall, there are a lot of remaining events around 
the edge of the fiducial volume. The reconstructed direc- 
tions of these remaining, externally-produced events in 
the fiducial volume are pointed, on average, strongly in- 
ward. We have eliminated most of these external events 
by using an event selection based upon vertex and direc- 
tion information. The details of this event selection will 
be explained in Section FVII Dl Although most of these 
external events are eliminated by this cut, some quantity 
of external events still remain. These remaining external 
events are one of the major backgrounds in this energy 
region. 

Another major source of background in the low energy 
region are radioactive daughter particles from the decay 
of ^^^Rn. ^^''Bi, which is one of these daughter particles, 
undergoes beta decay with resulting electron energies up 
to 3.26 MeV. Due to the limited energy resolution of the 
detector, these electrons can be observed in this energy 
region. 

We have reduced ^^^Rn in the water by our water pu- 
rification system to contribute less than 1 mBq/m'^ jl|, 
and monitor the radon concentration in real-time by sev- 
eral radon detectors [TolllTlIT^ . However, the water flow 
from the water inlets, located at the bottom of the SK 
detector, stirs radon emanated from the inner detector 
wall into the fiducial volume Therefore, there is an 
event excess after the external event cut in the detec- 
tor bottom region as compared to the top region. We 
also supply radon-reduced air into the space above 
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Isotope Ti (sec) Decay mode Kinetic Energy(MeV) 



^He 


0.119 


/9" 




9.67 + 0.98 ( 7 ) 








n 


( 16 % ) 


iLi 


0.838 


p~ 




~13 


IB 


0.77 


p~ 


h 


~13.9 


iLi 


0.178 


P' 




13.6 ( 50.5 % ) 






p~ 


n 


( -50 % ) 


Ic 


0.127 


p+ 


p 


3 ~ 15 




0.0085 


p~ 




16 ~ 20 (~50%) 








n 


-16 (-50%) 


i^Be 


13.8 


P~ 




11.51 ( 54.7 % ) 
9.41 + 2.1 ( 7 ) ( 31.4 % ) 


FBe 


0.0236 


P- 




11.71 




0.0202 


P~ 




13.37 




0.0110 


P" 


h 


16.32 


Pb 


0.0174 


P~ 




13.44 




0.0086 


P~ 


h 


13.2, 16.7 


^B 


0.0138 


P- 




14.55 + 6.09 ( 7 ) 


6 


2.449 


P- 




9.77 ( 36.8 % ) 
4.47 + 5.30 ( 7 ) 


16 

6 


0.747 


P- 


n 


~4 




7.13 


P- 




10.42 ( 28.0 % ) 
4.29 + 6.13 ( 7 ) ( 66.2 % ) 



TABLE V: Possible radioactive spallation products in Super- 
Kamiokande. 



the water surface in the SK tank to prevent radon in 
the mine air from dissolving into the purified water. The 
radon concentration in this radon-reduced air is less than 
3 mBq/m"^ and the measured radon level in the air in the 
tank is stable at 20~30 mBq/m^. Periods of anomalously 
high radon concentration due to water system troubles 
were removed from this analysis. 



B. High energy backgrounds 

Above about 6.5 MeV, the dominant background 
source is radioactive isotopes produced by cosmic ray 
muons' spallation process with oxygen nuclei. Some frac- 
tion of downward-going cosmic ray muons interact with 
oxygen nuclei in the water and produce various radioac- 
tive isotopes. These radioisotopes are called "spallation 
products." Table IVl shows a summary of possible spalla- 
tion products in SK. The (3 and/or 7 particles from the 
radioactive spallation products arc observed in Super-K, 
causing what are called "spallation events." The spal- 
lation events retain some correlation in time and space 
with their parent muons. Using this correlation, we have 
developed a cut to remove these spallation background 
events efficiently. The detail of this cut will be explained 
in Section lyilBl 



VII. DATA ANALYSIS 
A. Noise reduction 

The first step of the data reduction is an elimination 
of the noise and obvious background events. First of 
all, the events with total photoelectrons less than 1000, 
which corresponds to ~ 100 MeV, are selected. Next, 
the following data reductions are applied: (a) Events 
whose time difference since a previous event was less than 
50 /xsec were removed in order to eliminate decay elec- 
trons from stopping muons. (b) Events with an outer 
detector trigger and over 20 outer PMT hits were re- 
moved in order to eliminate entering events like cosmic 
ray muons. (c) A function to categorize noise events is 
defined by the ratio of the number of hit PMT's with 
IQI < 1.0 p.e. and the total number of hit PMT's. 
Since typical noise events should have many hits with 
low charge, events with this ratio larger than 0.4 were 
removed, (d) A function which can recognize an event 
where most of its hits are clustered in one electronics 
module [ATM] is defined. It is the ratio of the maximum 
number of hits in any one module to the total number 
of channels (usually 12) in one ATM module. Events 
with a larger ratio than 0.95 were removed as they gen- 
erally arose due to local RF noise in one ATM module, 
(e) Events produced by flasher PMT's must be removed. 
These flasher events often have a relatively larger charge 
than normal events, therefore they are recognized us- 
ing the maximum charge value and the number of hits 
around the maximum charge PMT. The slightly involved 
criteria is shown in Figure 1^ (f) An additional cut to 
remove noise and flasher events is applied using a com- 
bination of a tighter goodness cxii{goodness < 0.6) and 
a requirement of azimuthal uniformity in the Cherenkov 
ring pattern. Good events have a uniform azimuth dis- 
tribution of hit-PMT's along the reconstructed direction, 
while flasher events often have clusters of hit PMT which 
lead to a non-uniform azimuth distribution. 

In this noise reduction, the number of candidate solar 
neutrino events went from 3.43 x 10^ to 1.81 x 10^ after 
applying the 2 m fiducial volume cut and constraining 
the energy region from 5.0 to 20.0 MeV. Table IVIll shows 
the reduction step summary. 

The loss of solar neutrino signal by the noise reduction 
is evaluated to be 0.8% using a Ni(n,7)Ni gamma source 
at several edge positions of the inner detector. The differ- 
ence between data and MC for this source yields the sys- 
tematic error for this series of reductions; it's ±1.0% for 
our solar neutrino flux measurement. The loss is mainly 
due to the flasher cut. 



B. Spallation cut 

The method to reject spallation products shown in Sec- 
tion ^T] is described in this section. In order to identify 
spallation events, likelihood functions are defined based 
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FIG. 26: The relation between max charge and the number 
of hit PMT's in a 5 X 5 patch surrounding the PMT with the 
maximum charge. At most 24 tubes in this patch (plus the 
one at the center) can be activated. The upper plot shows 
a typical distribution including active flasher events, and the 
lower plot shows a typical good data set. The cut region is 
the area above the line in the upper figure. 



on the following parameters: 

• AL : Distance from a low energy event to the track 
of the preceding muon. 

• AT : Time difference from muon event to the low 
energy event. 

• Qres ■ Residual charge of the muon event, Q total — 
Qunit X Lfi, where Qtotai is the total charge, Qunit 
is the total charge per track length and is re- 
constructed track length of the muon track. 

Some fraction of muons deposit very large amounts of 
energy in the detector and in such cases vertex posi- 
tion reconstruction is not reliable. Therefore, the spalla- 
tion likelihood functions are defined for the following two 

cases, 

in the case of a successful muon track fit; 

Lspai^L, AT,Qres) = 

LfJ^,iAL,Qre.) X L%l{AT) X L%TiQres), (7.1) 

in the case of a failure fit the muon track; 

L^paiAT, Qtotai) = ifpI(AT) X L?-*- (Qtotai). (7.2) 

where Lfj,i{AL,Qres), ifpI(Ar), L%^riQres) and 
L'^pa'^^ [Qtotai) are likelihood functions for Ai, AT and 

Qres • 

Figure 1771 shows the AL distributions from spallation 
candidates for six Qres ranges, and the spallation-like 



function is made from these plots. Here, the selection 
criteria is AT < 0.1 sec and N^ff > 50 (equivalent to 
7.2 MeV). The peak around ~ 100 cm is caused by spal- 
lation events and that around 1500 cm is due to chance 
coincidence. The distribution of non-spallation events is 
calculated using a sample which shuffles the event times 
time randomly; it is the dashed line in the figure, and the 
non-spallation-like function is also made by these plots. 
After subtracting the non-spallation from the spallation 
function, and taking a ratio with a random coincidence 
distribution, Tf^^(AL, Q^es) was obtained. 
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FIG. 27: AL distribution for each Qrea ranges, (1) Qr 
2.4 X lO'* p.e. (2) 2.4 x 10'' < Q - ^ o in4 
4.8 X W < Qres < 9.7 X W p.e. 

4.8 X 10^ p.e. (5) 4.8 x 10^ < Qres < 9.7 x 10*^ p.e. (6) 
9.7 X 10'' p.e.< Qres- The solid line shows the data, and the 
dashed line shows the random sample. 



Figure EH shows the AT distributions from spallation 
candidates for each time range. Here, the selection crite- 
ria is AL < 300 cm and N^ff > 50 and Qres < 10^ p.e.. 
These distributions are fitted with the following function: 



-AT 
^spa 



(AT) 



7 

2^1 



AT 
-1/2 



(7.3) 



where t^^2 are half-life times of typical radioactive el- 
ements produced by spallation. The evaluated half-life 
times and radioisotopes are summarized in Table IVII 
In order to obtain the likelihood function for resid- 
ual charge L'^j^^" (Qres), time correlated events with low 
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FIG. 28: AT distribution for each time range, (1) 0.0 < AT < 
0.1 sec (2) 0.1 < AT < 0.8 sec (3) 0.8 < AT < 4.0 sec (4) 
4.0 < AT < 15 sec (5) 15 < AT < 100 sec (6) < AT < 
100 sec. Cross marks are data and lines are fitted likelihood 
function Lf^'^(AT) 




P-e- xlO 



1000 2000 3000^ 



FIG. 29: (a) Qres distribution for the spallation candidate 
events (solid line) and the non-candidate events (dashed line) . 
(b) Cross marks show the result of subtracting dashed from 
solid line in (a), and dotted line shows the likelihood function 

L^pa^ {Qres)- 



is calculated using a low energy sample which is shuffled 
in time and position. It has some position dependence 
because of the SK tank geometry. Figure 1^ shows the 
position dependence of the dead time as a function of the 
distance from the inner barrel wall and the top or bottom 
walls. The systematic error due to position dependence is 
estimated by comparing between MC and data, and it is 
estimated to be ±0.2% for flux, and ±0.1% for day-night 
flux difference and other time variations. 



i 


radioactivity 


''"1/2 


A, 


1 


FB 


2.02 X 10"^ 


1.20 X 10-' 


2 


FN 


1.10 X 10"^ 


3.39 X 10" 


3 


in 


1.78 X 10"' 


3.39 X 10^ 


4 


iLi 


8.40 X 10"' 


1.25 X 10^ 


5 


6 


2.45 


1.35 X 10^ 


6 


FN 


7.13 


6.76 X 10^ 


7 


FBe 


13.83 


7.79 



TABLE VI: The parameters of the hkelihood Tfpa(AT) 




Spallation likelihood 



Spallation likelihood 



energy events (AT < O.lsec, N^/f > 50) and non- 
correlated events are selected. Figure I29f a) shows the 
Qres distribution for spallation and non-spallation can- 
didates. After subtracting the non-spallation from the 
spallation function, the resulting distribution's fit by a 
polynomial function is shown in Figure IJHIb) . 

To employ the spallation cut, the likelihood values are 
calculated for all muons in the previous 100 seconds be- 
fore a low energy event, and a muon is selected which 
gives the maximum likelihood value(Lmaa;)- Figure QUI 
shows Lmax distributions for both the successful muon 
track reconstruction case and the failed reconstruction 
case. To be considered a spallation event the selec- 
tion criteria are Lmax > 0.98 (when fit succeeded) and 
Lmax > 0.92 (when fit failed). 

The dead time for low energy events caused by the 
spallation cut is estimated to be 21.1%. This estimation 



FIG. 30: The maximum likelihood value distribution for suc- 
cess (left) and failure (right) of muon track reconstruction. 
The blank histogram shows all data while the hatched his- 
togram shows the plot for the random sample; the latter is 
used to evaluate dead time. 



C. Ambient background reduction 

Even after the fiducial volume cut, some fraction of 
the ambient background still remains. It is mainly due 
to mis-reconstruction of the vertex position. In order to 
remove the remaining background, several cuts to evalu- 
ate the quality of the reconstructed vertex were applied. 

a. Fit stability cut The goodness value of the ver- 
tex was calculated for points in the region of the recon- 
structed vertex position and the sharpness of the good- 
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FIG. 31: Position dependence of the spallation cut's dead 
time. The horizontal axis shows the distance from the barrel 
(left), and the minimum distance from the top or bottom 
(right). 



ness distribution as a function of detector coordinates was 
evaluated. First of all, the goodness at ^ 300 grid points 
around the original vertex are calculated, and their differ- 
ence from the goodness at the reconstructed vertex (Ag) 
is determined. The number of test points which give a 
more than some threshold was counted; the threshold is 
as a function of energy and vertex position. The ratio of 
the number of points over this threshold to the total num- 
ber of grid points is defined as Rgrid- The Rgrid distri- 
bution of the data is shown in Figure 03 along with sim- 
ulated '^B Monte Carlo events. Events with Rgrid > 0.08 
are rejected as background. The systematic error of this 
reduction was evaluated using LINAC data and MC and 
the systematic error for flux measurement is estimated 
to be ±1.0%. 




FIG. 32: Rgrid distribution for data (solid line) and MC 
(dashed line). The normalization is done by scaling to the 
first bin. The dotted line shows the reduction criteria. 



b. Hit pattern cut Mis-reconstructed events often do 
not have the expected Cherenkov ring pattern when the 
hit PMT's are viewed from the reconstructed vertex. 
Also, some fraction of spallation products like ^^N emit 
multiple gammas in addition to an electron and so do not 



fit a single Cherenkov ring pattern very well. In order 
to remove those kind of events, the following likelihood 
function is defined: 



L{E x) EE E£illog(/(cos0rf,.,i?,f)). 



hit 



(7.4) 



here, f{cos9dirE, x) is the same function as that derived 
from MC in Eq. (14.31) , but here it depends on energy and 
vertex position. Figure 123 shows the likelihood distribu- 
tion Eq. (fn|l for data and MC. The events with a 
likelihood less than —1.85 are rejected. The systematic 
error of this cut is estimated using solar neutrino MC 
and a sample of very short lived spallation products; it 
is -1-1.0% and —0.5% for the flux measurement. 
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FIG. 33: Likelihood distribution for data (solid line) and MC 
(dashed line). The MC is normalized to the peak bin. The 
dotted line shows the reduction criteria. 



c. Fiducial volume cut using the second vertex fitter 
The 2 m fiducial volume cut using the second vertex re- 
construction described in Section IIV Al is applied. Fig- 
ure 1^ shows the distance from the wall distribution of 
data and MC using the second vertex reconstruction. Ta- 
ble I VIII also shows the reduction step summary described 
in this section. 



D. Gamma ray cut 

The gamma-ray cut is used for reduction of external 
gamma rays which mainly come from the surrounding 
rock, PMT glass, and stainless steel support structure of 
the detector. These gamma rays are one of the major 
background especially for solar neutrino data. 

The distinctive feature of external gamma rays is that 
they travel from the outside edge of the SK volume and 
continue on to the inside. In order to remove this kind of 
event, reconstructed direction and vertex information is 
used to define the effective distance from the wall, de//, 
as shown in the inset of the left plot in Figure As 



19 



X 10 ' 




50 100 150 200 250 300 350 400 
distance from wall (cm) 



FIG. 34: The distance from the wall distribution using 
the second vertex reconstruction for data (solid) and MC 
(dashed) . The normalization of MC is done by the total num- 
ber of events. The dotted line shows the reduction criteria. 



they emanate from the wall itself, the value de// for ex- 
ternal gamma ray events should be small. Figure IHBl also 
shows de// distributions for data and MC. The criterion 
of the gamma ray cut is determined by maximizing its 
significance; here the number of remaining events after 
reduction in data and MC are used to define the signif- 
icance, NMc/v^data- The gamma ray cut criteria are 
determined to be: 

1. de// > 450 cm {ioT E> 6.5 MeV ) 

2. de// > 800 cm ( for 5.0MeV < £; < 6.5 MeV) 



Since the gamma ray cut uses reconstructed vertex and 
direction, the differences of vertex and angular resolution 
between data and MC can introduce systematic errors. 
In order to estimate this effect, we shift the reconstructed 
vertex and direction of events within the difference of 
data and MC, which are measured by LIN AC, and apply 
the gamma ray cut to these events. The difference of 
reduction efficiency before and after this shift is taken 
to be the systematic error for the gamma ray cut. The 
systematic error of this cut for the flux measurement is 
estimated to be ±0.5% , while for the energy spectrum 
it is estimated to be ±0.1%. 
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FIG. 35: Effective distance (d^//) oi E > 6.5 MeV (left) and 
5.0MeV < E <6.5 MeV (right). Blank histogram shows data 
and hatched area shows solar neutrino MC. 



Figure and |27| show the vertex and direction dis- 
tributions before and after the gamma ray cut for real 
data in the energy regions E > 6.5 MeV and 5.0 MeV 
< E < 6.5 MeV. The dead times introduced by this cut 
in the energy regions i? > 6.5 MeV and 5.0 MeV < E < 
6.5 MeV are estimated to be 6.9% and 22.0%. 



FIG. 36: Vertex distribution before (blank) and after 
(hatched) applying the gamma ray cut. The left plots are 
for |Z| <1000 cm, and the right ones are Z for |R| <1000 cm. 
The upper plots show E > 6.5 MeV, and lower ones show 
5.0 MeV < E <6.5 MeV. 



^^N background is generated by the capture of /i on 
^^O nucleus in water: 

The most probable decay mode of produces a 

6.1 MeV 7 ray together with a (3 decay electron of max- 
imum energy 4.3 MeV; its half life is 7.13 seconds. 

To tag this kind of event, spatial and time correlations 
with captured stopping muons are used. The selection 
method is as follows: 1) in order to select only a captured 
muon sample, collect a sample of stopping muons which 
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FIG. 37: Directional distribution (X, Y, Z) before (blank) and 
after (hatched) applying the gamma ray cut. The upper plots 
show E > 6.5 MeV, and the lower plots show 5.0 MeV < E < 
6.5 MeV. 
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FIG. 38: Summary of the reduction efficiencies on MC. 



are not followed by a decay event in 100 usee, then 2) 
select a low energy event sample within 335 cm from the 
stopping point of the muon as well as in a time window 
of 100 msec to 30 seconds following the stopping muon. 
The number of such candidate low energy events is 9843 
in 1496 days of low energy data. The dead time for this 
reduction is estimated by using a so-called random sam- 
ple, where the time and vertex positions of the events 
have been mixed randomly. The same ^^N event selec- 
tion as described above is then applied to the random 
sample, yielding a dead time for the ^^N cut of 0.49%. 



F. Event reduction summary 

Figure |2H1 shows the reduction efficiencies after each 
step as a function of energy using MC simulation events. 
The efficiency is used as a correction when the energy 
spectrum is calculated as described in Section IVIII Dl 
Figure 123 shows the event rate after each reduction step 
for data and also shows the predicted solar neutrino spec- 
trum. Table IVni summarizes the numerical results of the 
reduction steps. The number of events after all the re- 
duction steps is 286557. 



VIII. RESULTS 

A. Solar neutrino signal extraction 

The solar neutrino signal is extracted from the direc- 
tional correlation of the recoiling electrons with the inci- 
dent neutrinos in v-e scattering. FigureHUIshows cos ^sun 
where ^sun is the angle between the reconstructed re- 
coil electron direction and the expected neutrino direc- 
tion (calculated from the position of the sun at the event 
time). 
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FIG. 39: Summary of the data reduction steps. 



The solar neutrino flux is extracted by a likelihood 
fit of the solar neutrino signal and the background to 
this distribution. This likelihood function is defined as 
follows: 

C = e~(^' ^'+^) n n • + S-Y^- (8.1) 

i=l 3 = 1 

We define A'^tin = 21 energy bins: 18 energy bins of 
0.5 MeV between 5.0 and 14.0 MeV, two energy bins 
of 1 MeV between 14.0 and 16.0 MeV, and one bin be- 
tween 16.0 and 20.0 MeV. S is the total number of solar 
neutrino signal events, and n^, Bi, and Yi represent the 
number of observed events, the number of background 
events, and the expected fraction of signal events in the 
i-th energy bin, respectively. We use two types of proba- 
bility density functions: p(cos6'aun, E) describes the angu- 
lar shape expected for solar 's of recoil electron energy 
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Total 


3.43 X 10' 


A. 


Noise reduction 


(a) 


2.66 X 10^ 


(b) 


2.51 X 10^ 


(c) 


2.50 X 10^ 


(d) 


2.50 X lO'' 


(e) 


2.48 X 10'' 


(f) 


1.81 X 10^ 


B. 


Spallation cut 




1.29 X 10^ 


C. 


Ambient B.G. cut 


(a) 


3.61 X 10*^ 


(b) 


2.72 X 10® 


(c) 


1.86 X 10*^ 


D. 


Gamma cut 




2.96 X 10^ 


E. 


^"N cut 




2.87 X 10^ 



TABLE VII: The summary of number of events remaining 
after each reduction step 



E (signal events) and (cos 6*500) is the background shape 
in energy bin i. Each of the rii events in energy bin i is 
assigned the background factor bij = m (cos 9ij ) and the 
signal factor Sij = p{cos9ij, Ej). 

The signal shape p(cos6'sun, i?) is obtained from the 
known, strongly forward-peaked angular distribution of 
neutrino-electron elastic scattering with smearing due to 
multiple scattering and the detector's angular resolution. 
The background shape Ui(cos0sun) has no directional cor- 
relation with the neutrino direction, but deviates from a 
flat shape due to the cylindrical shape of the SK de- 
tector: the number of PMT's per solid angle depends 
on the SK zenith angle. In order to calculate the ex- 
pected background shape, we use the angular distribu- 
tion of data itself. The presence of solar neutrinos in the 
sample biases mostly the azimuthal distribution, so at 
first we fit only the zenith angle distribution and assume 
the azimuthal distribution to be flat. We generate toy 
Monte Carlo directions according to this fit and calcu- 
late cos^Bun. We also fit both zenith and azimuthal dis- 
tributions, approximately subtracting the solar neutrino 
events from the sample and repeat the toy Monte Carlo 
calculation. We compare the obtained number of solar 
neutrino events from both background shapes and as- 
sign the difference as a systematic uncertainty. Since the 
azimuthal distributions don't deviate very significantly 
from fiat distributions, we quote the solar neutrino events 
obtained from the first shape (assuming a fiat azimuthal 
distribution) . The dotted area in Figure QUI shows this 
background shape. The systematic uncertainty due to 
the background shape is 0.1% for the entire data sample 
(5.0-20.0 MeV). If the data sample is divided into a day 
and a night sample, the systematic uncertainty is 0.4%. 
The amount of background contamination is much less 
above 10 MeV than it is near the SK-I energy thresh- 
old (5.0 MeV), so small differences in background shape 



between the two methods become important only in the 
lowest energy bins: between 5.0 and 5.5 MeV, the sys- 
tematic uncertainty is estimated to be 1.2%, between 5.5 
and 6.0 MeV 0.4%, and above 6.0 MeV 0.15%. 




FIG. 40: Angular distribution of solar neutrino event candi- 
dates. The shaded area indicates the elastic scattering peak. 
The dotted area is the contribution from background events. 



B. Observed solar neutrino flux 

Figure shows the cos ^sun distribution for 1496 days 
of SK-I data. The best fit value for the number of 
signal events due to solar neutrinos between 5.0 MeV 
and 20.0 MeV is calculated by the maximum likeli- 
hood method in Eq. (|8.1() . and the result for SK-I is 
22, 404 ± 226 (stat.)+^^^ (sys.). The corresponding 
flux is: 

2.35 ± 0.02 (stat.) ± 0.08 (sys.) x 10*^ cm'^s'^ 

The systematic errors for the solar neutrino flux, sea- 
sonal variation and day-night differences for the energy 
range 5.0 MeV to 20.0 MeV are shown in Table lvml The 
detailed explanations are written in each topic's section, 
but the total systematic error for the solar neutrino flux 
measurement is estimated to be ^o ^S. 



C. Time variations of solar neutrino flux 

1. Day-Night difference 

The day time flux and night time flux of solar neutrinos 
in SK-I are calculated using events which occurred when 
the solar zenith angle cosine was less than and greater 
than zero, respectively. The observed flux are: 

^day = 2.32 ±0.03 (stat.)j:°:°^ (sys.) x 10^ cm-^s"! 
<^mght = 2.37 ±0.03 (stat.)t°:°^ (sys.) x 10^ cm^^^-i 
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Flux Seasonal day-night 


Energy scale, resolution 


±1.6 


+ 1.2 

— 1.1 


+ 1.2 
— 1.1 


Theoretical uncertainty 


+1.1 

— 1.0 






for *B spectrum 








Trigger efficiency 


+0.4 
—0.3 


±0.1 




Reduction 


+2.1 
— 1.6 


±0.5 




Spallation dead time 


±0.2 


±0.1 


±0.1 


Gamma ray cut 


±0.5 


±0.25 




Vertex shift 


±1.3 






Background shape 


±0.1 




±0.4 


for signal extraction 








Angular resolution 


±1.2 






Cross section of I'-e scattering 


±0.5 






Livetime calculation 


±0.1 


±0.1 


±0.1 


Total 


+3.b 
-3.2 


±0.3 


+ 1.3 
-1.2 



TABLE VIII: Systematic error of each item (in %) 



Their difference leads to a day-night asymmetry, defined 

as ^ = {^day - ^7itght)/{^{'^day + '^mght))- We find: 




A 



-0.021 ±0.020 (stat.)l[;;^}^ (sys.) 



FIG. 41: The solar zenith angle dependence of the solar neu- 
trino flux (error bars show statistical error). The width of 
the night-time bins was chosen to separate solar neutrinos 
that pass through the Earth's dense core (the rightmost Night 
bin) from those that pass through the mantle. The horizontal 
line shows the flux for all data. 



Including systematic errors, this is less than 1 — cr from 
zero asymmetry. The largest sources of systematic error 
in the asymmetry are energy scale and resolution (Ig on) 
and the non-flat background shape of the cos Osun distri- 
bution (±0.004). As described in the neutrino oscillation 
analysis section, we can reduce the statistical uncertainty 
if we assume two-neutrino oscillations within the Large 
Mixing Angle region. The day-night asymmetry in that 
case is 

A = -0.017 ± 0.016 (stat.)t[!;!]}^ (sys.) ± 0.0004 (osc.) 

with the final, tiny additional uncertainty due to the un- 
certainty of the oscillation parameters themselves. Fig- 
ure ^2 shows the solar neutrino flux as a function of the 
solar zenith angle cosine. 
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2. Seasonal variation 

Figure 021 shows the monthly variation of the flux, 
which each horizontal bin covers 1.5 months. The flg- 
urc shows that the experimental operation is very stable. 

Figure 021 shows the seasonal variation of solar neu- 
trino flux. As in Figure 021 each horizontal time bin is 
1.5 months wide, but in this figure data taken at simi- 
lar times during the year over the entire course of SK-Fs 
data taking has been combined into single bins. The 1.7% 
orbital eccentricity of the Earth, which causes about a 
7% flux variation simply due to the inverse square law, 
is included in the flux prediction (solid line). The ob- 
served flux variation is consistent with the predicted an- 
nual modulation. Its x^/d.o.f. is 4.7/7, which is equiva- 
lent to 69% C.L.. If we fit the eccentricity to the Earth's 
orbit to the observed SK rate variation, the perihelion 



FIG. 42: Solar neutrino flux as a function of time. The bin- 
ning of the horizontal axis is 1.5 months. 



shift is 13 ± 18 days (with respect to the true perihelion) 
and the eccentricity is 2.1 ±0.3% 14]. This is the world's 
first observation of the eccentricity of the Earth's orbit 
made with neutrinos. The total systematic error on the 
relative flux values in each seasonal bin is estimated to 
±1.3%. The largest sources come from energy scale and 
resolution (i*^}'^^) and reduction cut efficiency (±0.5%), 
as shown in Table lyTTll 



D. Energy spectrum 

Figure^lshows the expected and measured recoil elec- 
tron energy spectrum. The expected spectrum is calcu- 
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FIG. 43: The seasonal variation of the solar neutrino flux. 
The solid line is the prediction based on the eccentricity of 
the Earth's orbit. 



lated by the detector simulation described in Section HIH 
and BP2004 is used as a solar model. The solid line shows 
the expected spectrum of *B and hep neutrinos, and the 
dashed line shows the contribution of only neutrinos. 

The observed and expected event rates are summarized 
in Table ITXI and Table 0C| In these tables the reduction 
efficiencies listed in Fig. |2H1 are corrected. 

The uncorrelated and correlated systematic errors for 
each energy bin are shown in Table EH and Table IXTTI 
respectively. 
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FIG. 44: Energy spectrum of the solar neutrino signal. The 
horizontal axis is the total energy of the recoil electrons. The 
vertical axis is the event rate of the observed solar neutrino 
signal events. The error bars are a quadrature of the sta- 
tistical and uncorrelated errors. The reduction efficiencies in 
Fig. EH are corrected. BP2004 SSM flux values are used for 
the *B and hep MC fluxes in this plot. The dashed line shows 
the contribution of only *B. 



"10 
(D ' ^ p 

LU10 

1-1 
10 

■5 3 

cc 

10^ 



10 





Integrated up to 21 MeV ■ 


r 


B8+HEP 90%UL : 




only hep' I ■ 


1 I 1 1 1 


1 1 1 1 1 1 1 I 1 1 1 1 1 1 1 1 1 1 1 I 1 
— 1 Best aperture 




' — 1 = 18-21 MeV ^ 



8 10 12 14 16 18 20 
Energy thresold(MeV) 

FIG. 45: Integral energy distributions for the hep solar neu- 
trino analysis. The horizontal axis is the energy threshold of 
the recoil electrons. The integration is up to 21.0 MeV. Up- 
per plot: the vertical axis is the expected number of events 
in SK-I's 1496 day final data sample after all cuts. The solid 
and dashed lines correspond to the hep solar neutrino only 
the and *B-|-hep 90% C.L. upper limit. Lower plot: the ratio 
of the two histograms in the upper plot (dashed/solid). 



E. Hep solar neutrino 



The expected hep neutrino flux is three orders of mag- 
nitude smaller than the ^B solar neutrino flux. However, 
since the end-point of the hep neutrino spectrum is about 
18.8 MeV compared to about 16 MeV for the ^B neutri- 
nos, the high energy end of the *B spectrum should be 
relatively enriched with hep neutrinos. In order to dis- 
cuss the flux of hep neutrinos, the most sensitive energy 
region for hep neutrino was estimated, then, assuming 
all the signal events in this energy region were due to 
hep neutrinos, an upper limit of the hep solar neutrino 
flux was obtained. Any possible effects from neutrino 
oscillations were not considered in this analysis. 

Figure shows the expected integral energy distribu- 
tions for solar neutrinos. In the high energy region, the 
relative hep contribution is high, but the expected num- 
ber of events is small because of the limited observation 
time. For this analysis, the best energy aperture for the 
hep solar neutrinos was determined to be 18.0~21.0 MeV. 
In this energy region 0.90 hep solar neutrino events are 
expected from the predicted BP2004 SSM rate. 

Applying the same signal extraction method to the 
data events in the 18.0'--^21.0 MeV region, we found 
4.9 ± 2.7 solar neutrino signal events. Assuming that 
all of these events are due to hep neutrinos, the 90% 
confldence level upper limit of the hep neutrino flux was 
73 X 10^ cm~^ s~^. Figure 021 shows the differential en- 
ergy spectrum of solar neutrino signals in the high energy 
region. 
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Energy 




Observed rate 




Expected rate 


(MeV) 


ALL 


DAY 


NIGHT 


B 


hep 






-1 < COS^z < 1 


-1 < cos 6*2 < 


< cos 6iz < 1 






5.0 


-5.5 


74.7^6:? 


72.1^9;:^ 


77 i+y.2 

' '-^^-9.0 


182.9 


0.312 


5.5 


-6.0 


65.0t|J 


64.811^ 


65.ltlj 


167.7 


0.309 


6.0 


- 6.5 


61.5t^j 


60.211^ 


62.6^^3 


151.9 


0.294 


6.5 


- 7.0 


54.lti:? 


54.21^4 


53.9t^j 


135.3 


0.284 


7.0 


- 7.5 


49.4+^-^ 


49.2+1? 


49.6t|.l 


119.2 


0.266 


7.5 


-8.0 


44.3t;i 


44.8+?'° 


A3.8t\i 


103.5 


0.249 


8.0 


-8.5 


36.3+^-2 


35.7lp7 


36.8^^7 


88.3 


0.236 


8.5 


-9.0 


28.7+i-« 




30.6^^-4 


74.1 


0.221 


9.0 


-9.5 




25.4tl-j 


24.6tiJ 


61.4 


0.196 


9.5 


- 10.0 


20.8t°:« 


2O.7+I1 


20.8+ii 


49.9 


0.185 


10.0 


- 10.5 


16.2+0-^ 


1 5 7+1'° 
-0.9 


16 7+^'-° 

' -0.9 


39.6 


0.167 


10.5 


- 11.0 


1 1 r,+0.6 


10.9+0? 


11-O-0.8 


30.7 


0.149 


11.0 


- 11.5 


9.85_o;49 


65+0 '''* 


10.031°;^! 


23.28 


0.130 


11.5 


- 12.0 


6 70+0-42 

' •'-0.40 


7 14+0.63 


6 47+0-58 


17.27 


0.118 


12.0 


- 12.5 


c 1 n+0.36 
^^'^-0.33 


c 91+0.51 
"-••^^-0.46 


12.45 


0.098 


12.5 


- 13.0 


o f;r+0.30 


3.96tr4? 


q qo+0.41 
0.00_o.3g 


8.76 


0.090 


13.0 


- 13.5 


^■*O-0.23 


2.56«:^3« 


9 07+0.35 
^■'^'-0.30 


5.94 


0.073 


13.5 


- 14.0 


o rio + 0.22 
^•"^-0.20 


1 95+0-32 
i.yo_g 27 


9 00+0-^2 
z.u»_(j 27 


3.88 


0.060 


14.0 


- 15.0 


1 79+0.21 
J^- '^-0.19 


l.DU_o 25 


1 fic;+0-31 


4.01 


0.094 


15.0 


- 16.0 




1 1 q/^+0.238 

1.10U_Q j^g2 


1.439 


0.057 


16.0 


- 20.0 


Q4I +0.103 


u.zy:U_(j 093 


49'^+0-156 


0.611 


0.068 



TABLE IX: Observed and expected event rates in each energy bin at 1 AU. The unit of the rates is events/kton/year. The 
errors in the observed rates are statistical only. The reduction efficiencies in Fig. l38l are corrected, and the expected event rates 
are for the BP2004 SSM flux values (*B flux is 5.79 xlO'', hep flux is 7.88 x 10^/cm^/sec). 6^ is the angle between the z-axis 
of the detector and the vector from the Sun to the detector. 



Energy 
(MeV) 


DAY 

cos 6*2 : — 1 ~ 


MANTLEl 
~ 0.16 


Obf 

MANTLE2 
0.16 ~ 0.33 


erved rate 
MANTLE3 
0.33 ~ 0.50 


MANTLE4 
0.50 ~ 0.67 


MANTLES 
0.67 ~ 0.84 


CORE 
0.84 ~ 1 


Expec 
»B 


;ted rate 
hep 


5.5-6.5 
6.5 - 8.0 
8.0-9.5 
9.5 - 11.5 
11.5 - 13.5 
13.5 - 16.0 


127.+^- 
149.tl 

57.lt\i 
I8.7tj;° 
4.28l°:« 


124.:-; 
166.t-: 
90.7tli 

20.0+24 

4.78ll;« 


106.til- 
158.+f- 
92.1+11 
63.3ti° 
13.8t^-.g 

6.97t\il 


132.^^2. 
137.tg' 
90.5+^-^ 
56.8+^9 
15.3t?;° 

p: 09+1.22 
•^•"^-0.98 


146.ii^; 
150.tr. 
99.8t^:^g 
59.61^9 
19.5tll 

r 50 + 1.19 


140. +i^' 

141. +«g' 
90.3t^:^ 

60.ltt3 
17.0t^;^ 
70+1.14 


iw.lit 

137.+f- 
88.5t^;° 
60.9tl^ 
20.4t^-.^ 

Q qo + 1.27 

0.»0_(j g3 


320. 
358. 
223.8 
143.4 
44.4 
9.33 


0.603 
0.799 
0.653 
0.631 
0.379 
0.211 



TABLE X: Observed and expected event rates in each zenith-spectra data set at a one astronomical unit [AU] distance from the 
Sun. The unit of the rates is events/kton/year. The errors in the observed rates are statistical only. The reduction efficiencies 
in Fig. 1381 are corrected, and the expected event rates are for the BP2004 SSM flux values. 9^ is the angle between z-axis of 
the detector and the vector from the Sun to the detector. 



IX. SOLAR NEUTRINO OSCILLATION 
ANALYSIS 

A. Introduction 

Two-neutrino oscillations are so far sufficient to ex- 
plain and describe all measured solar neutrino phenom- 
ena. The flavor eigenstates Ve and (where is either 

or i^r or an admixture of both) describe weak interac- 
tions of neutrinos and electrons or nucleons. Only solar 
j/g's can participate in charged-current reactions, since 
the solar neutrino energy is below the mass threshold. 



These flavor eigenstates are related to the mass eigen- 
states vi and 1^2 via the unitary mixing matrix U, which 
for two neutrinos can be expressed in terms of a single 
parameter, the weak mixing angle 6: 

jj _ f cos 6 sin 6\ 
\ — sin 9 cos 9 J 

In vacuum, the neutrino wavefunction oscillates in space 
with a frequency of \/p'^ + ~ p-\-m? /2p leading to an 
oscillatory transition probability of the flavor eigenstates 
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Energy (MeV) 


5 - 5.5 


5.5-6 


6-6.5 


6.5 - 7 


7- 


Trigger efficiency 


+'2.b 
-1.5 


-0.6 


±0.1 


±0.2 




Reduction 


±0.9 


±0.9 


±0.9 


±0.9 


±0.9 


Gamma ray cut 


±0.1 


±0.1 


±0.1 


±0.1 


±0.1 


Vertex stiift 


±0.2 


±0.2 


±0.2 


±0.2 


±0.2 


Background stiape 


±0.6 


±0.5 


±0.1 


±0.1 


±0.1 


for signai extraction 












Anguiar resoiution 


±2.3 


±1.0 


±1.0 


±1.0 


±1.0 


Cross section of 


±0.2 


±0.2 


±0.2 


±0.2 


±0.2 


iz-e scattering 












Total 


+3.b 
-29 


+ 1.V 

-1.6 


±1.4 


±1.5 


±1.4 



TABLE XI: Energy uncorrelated systematic errors for eadi 
energy bin. 



Energy 


Scale(%) 


Resolution(%) 


Tfieory(%) 


5.0 


- 5.5 


±0.1 


0.0 


±0.2 


-0.2 


±0.1 


0.0 


5.5 


-6.0 


-0.1 


±0.1 


±0.2 


-0.2 


0.0 


±0.1 


6.0 


-6.5 


-0.3 


±0.2 


±0.2 


-0.2 


-0.1 


±0.1 


6.5 


- 7.0 


-0.5 


±0.4 


±0.2 


-0.2 


-0.3 


±0.2 


7.0 


- 7.5 


-0.7 


±0.6 


±0.2 


-0.2 


-0.4 


±0.4 


7.5 


-8.0 


-0.9 


±0.9 


±0.2 


-0.2 


-0.5 


±0.5 


8.0 


-8.5 


-1.1 


±1.1 


±0.2 


-0.2 


-0.7 


±0.7 


8.5 


-9.0 


-1.4 


±1.4 


±0.1 


-0.1 


-0.9 


±0.9 


9.0 


-9.5 


-1.7 


±1.8 


±0.1 


-0.1 


-1.1 


±1.1 


9.5 


- 10.0 


-1.9 


±2.1 


0.0 


0.0 


-1.3 


±1.4 


10.0 


- 10.5 


-2.3 


±2.5 


-0.1 


±0.1 


-1.5 


±1.7 


10.5 


- 11.0 


-2.6 


±2.8 


-0.3 


±0.2 


-1.8 


±2.0 


11.0 


- 11.5 


-3.0 


±3.2 


-0.5 


±0.4 


-2.1 


±2.3 


11.5 


- 12.0 


-3.4 


±3.6 


-0.8 


±0.7 


-2.4 


±2.6 


12.0 


- 12.5 


-3.8 


±4.1 


-1.1 


±0.9 


-2.7 


±3.0 


12.5 


- 13.0 


-4.3 


±4.5 


-1.4 


±1.3 


-3.1 


±3.3 


13.0 


- 13.5 


-4.8 


±5.0 


-1.9 


±1.7 


-3.4 


±3.8 


13.5 


- 14.0 


-5.4 


±5.5 


-2.4 


±2.2 


-3.8 


±4.2 


14.0 


- 15.0 


-6.3 


±6.3 


-3.3 


±3.1 


-4.4 


±5.0 


15.0 


- 16.0 


-7.7 


±7.6 


-4.9 


±4.7 


-5.1 


±6.2 


16.0 


- 20.0 


-9.9 ±10.2 


-8.1 


±8.2 


-5.6 


±8.7 



TABLE XIL Energy correlated systematic errors for each of 
the 21 energy bins. 



>,10 

CO 

iio 

w 

|io 

(D 

T3 
(D 

lo 
o 

10 



1 



°B: 2.35x10 + hep: 73x10 
''6:2.35x10'^ + hep: 



10 



12 



14 



16 



18 20 
Energy (MeV) 



FIG. 46: Energy spectrum of recoil electrons in the high en- 
ergy region. The points show data with statistical error bars. 
The curves show expected spectra with various hep contri- 
butions to the best-fit ®B spectrum. The unit of the fluxes 
for the curves are /cm^/s. The dotted and dashed curves 
show the spectrum with 73 and (xlO^/cm^/s) hep fluxes, 
respectively. 
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In a 



where the oscillation length iosc = 
medium of matter density p, the small angle scattering 
of the Ve flavor differs (due to the additional charged- 
current amplitude) from the flavor; this can be de- 
scribed with a matter potential AV{p). The transition 
probability is still given by Eq. (|9.1|) but with the "effec- 
tive" oscillation length and mixing angle 



■^eff = 



TT 



1.27A/(2Ay - 



^ cos 26 



(^ sin 26*)^ 



FIG. 47: The Sun's electron density profile from the stan- 
dard solar model m [ssM] as a function of the distance to 
the center (top panel), and Earth's mass density profile (bot- 
tom panels). The left plot shows the density as a function of 
the distance to the Earth's center. The right plot shows the 
minimum (dashed horizontal line), average (solid line) and 
maximum density (dashed line) a solar neutrino of zenith an- 
gle 6z encounters on its path through the Earth. The electron 
density is obtained from this by multiplying by 0.497 mol/cm^ 
for the mantle and 0.468 mol/cm'^ for the core. 



tan 20gff 



■ sin 29 



2AV - 



^ cos 261 



These relations are valid only for a constant matter po- 
tential (and therefore constant matter density). We cal- 
culate the oscillation probability by means of a numerical 
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simulation using the position-dependent matter density 
profile in the Sun and Earth (shown in Figure 1^71) . 

As noted by Mikheyev, Smirnov, and Wolfenstein 
[MSW] [III, the matter density in the Sun (see Figure lT7|) 

is large enough for a resonance {AV — — cos 20) to 
occur, producing effective maximal mixing (6'gg- — j) 
even if the fundamental mixing 9 is small. The so- 
lar and terrestrial matter densities influence the tran- 
sition probability for solar neutrinos between Am^ /E w 
10-9 eVVMeV and « 10"^ eVVMeV [MSW range]. For 
the so-called Large Mixing Angle [LMA; 6* « 7r/6] MSW 
"solution" to the solar neutrino problem, Arri^ is chosen 
so that the solar neutrino spectrum lies at the "end" of 
this range (near 10"^ eV^/MeV). It explains all observed 
solar neutrino interaction rates. In the case of the Small 
Mixing Angle [SMA; 9 « 0.05] MSW solution, which also 
explains all observed rates, the solar neutrino spectrum 
must be placed near lO""^ eV^/MeV (closer to the center 
of the MSW range). Other large angle solutions exist as 
well: the Low Am^ solution [LOW] lies near the "be- 
ginning" of the MSW range (10"^ eV^/MeV), while the 
(quasi-) Vacuum [VAC] solutions are "below" the MSW 
range (lO-^^ eVVMeV to 10"" eVVMeV). 

To calculate the solar neutrino interaction rate on 
Earth, three steps are required: (i) the probability pi 
{P2) of a solar neutrino, which is born as Ve in the core of 
the Sun, to emerge at the surface as vi {^2), (ii) coherent 
or incoherent propagation of the vi , 1^2 admixture to the 
surface of the Earth, (iii) the probability pie {p2e) for a vi 
(1^2) neutrino to appear as Ve in the detector (after prop- 
agation through part of the Earth, if the Sun is below 
the horizon, using the PREM density profile of the 
Earth as shown in Figure ITTjl . In and above the MSW 
range of Arn^/E, the distance between the Sun and the 
Earth is much larger than the vacuum oscillation length, 
so the propagation (ii) can be assumed to be incoherent. 
In that case, the total survival probability of the Ve flavor 
is 



Pe 



■ Pi X Pie +P2^ P2e = "^PlPle + I - Pi - Pie 



where pi and pie are computed numerically. Below the 
MSW range, both the solar and terrestrial matter densi- 
ties can be neglected. However, the distance L between 
the Sun and the Earth approaches the oscillation length, 
so (ii) must be done coherently, and the survival proba- 
bility of the flavor is 



Pe 



1 -Pe 



1 - sin^ 29 81-0^ ( TT- 



L 



using Eq. (|9.1|l . Figure B51 shows (as an example) the 
survival probability of *B neutrinos. Other solar neu- 
trino branches may have slightly different probabilities, 
because the radial distributions of the neutrino produc- 
tion location differ. 

Neutrino oscillations impact SK data in three inde- 
pendent ways: (i) Since electron neutrinos have a much 
larger elastic scattering cross section than other flavors. 
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FIG. 48: Survival probability for B solar neutrinos. The top 
plot displays contours of equal probability as a function of 
Am^ /E and tan^ 6. The middle panel shows the same proba- 
bility as a function of AE/m^ only where a particular (small) 
mixing angle has been selected. The shaded area contrasts 
the difference between neutrinos passing through the core of 
the Earth (which has the highest matter density) with neutri- 
nos arriving directly from the Sun. The bottom panel shows a 
large mixing case. The line above the shaded area depicts the 
average probability for neutrinos passing through the Earth. 
Superimposed on these two lower panels are the locations of 
three neutrino branches (pp - dotted, ^Be - dashed, and *B 
- solid) for the SMA, LMA, LOW, VAC, and quasi-VAC so- 
lutions. 



neutrino oscillations reduce the rate of solar neutrino 
interactions, (ii) The spectrum of recoil electrons is 
distorted due to the energy dependence of the survival 
probability, (iii) The influence of the Earth's matter on 
the survival probability induces an apparent time depen- 
dence of the solar neutrino interaction rate with a 24 
hour period. The amplitude of that time dependence is 
expressed in the day/night asymmetry q c?n"i^A7^ where 



i.5{D+N] 

D (N) is the averaged interaction rate during the day 
(night). Due to the eccentricity of the Earth's orbit the 
distance between Sun and Earth changes periodically. 
Since the survival probability depends on this distance 
(if the neutrinos propagate coherently), this leads to an- 
other time- variation (with a 365 day period) expressed in 
the summer/winter asymmetry q 5(s+w) '^^^re S (W) is 
the averaged interaction rate during the summer (win- 
ter), corrected for the l/r^ dependence of the solar neu- 
trino intensity. The day/night and summer/winter varia- 
tions cannot both be present in the SK data, since above 
a Am^ of about 10^® eV'^ neutrinos propagate incoher- 
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FIG. 49: Expected SK day/night asymmetry (left) and spec- 
tral distortion (right) in the MSW region. The scale of the 
contours of equal asymmetry (distortion) is logarithmic. The 
spectral distortion is measured in terms of a ^ using the SK 
spectrum uncertainties (but not the values of the SK spec- 
trum). The gray-shaded areas correspond to the LMA, SMA, 
and LOW solutions. 



FIG. 51: Distortion of the recoil electron spectrum using 21 
(upper panel) and 8 (lower panel) bins. The measured event 
rates in each bin are divided by the event rates expected from 
MC assuming a *B flux of 5.076 x 10^/cm^s and a hep v flux 
of 32x 10^/cm'^s. Overlaid is a typical LMA solution (red) and 
a LOW solution (green), where the assumed neutrino fluxes 
were fit to best describe the data. 
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FIG. 50: Expected SK summer/winter asymmetry (left) and 
spectral distortion (right) in the vacuum oscillation region. 
The scale of the contours of equal asymmetry (distortion) is 
logarithmic. The spectral distortion is measured in terms of 
a using the SK spectrum uncertainties (but not the values 
of the SK spectrum). The gray-shaded areas correspond to 
the VAC solutions. 



ently and below 10~® eV^ the Earth's matter effects on 
the survival probabihty are negligible. Figure depicts 
the expected SK day /night asymmetry and FigurelCTthe 
summer/winter asymmetry, depending on the oscillation 
parameters. They also show the SK sensitivity to distor- 
tions of the recoil electron spectrum. 

Using the survival probability Pe, the neutrino inter- 
action rate due to the neutrino spectrum l(Ev) at SK 



is 

ro.c= \ dE I dE,I{E,) f dE,R{Ee,E){Se{E,,E,)p,+ 

S^{E,,E,){\-p,)) 

where Se.x{Eu, Eg) describe the probability that the elas- 
tic scattering of a Ve,x of energy E^ with electrons 
produces a recoil electron of energy Ee- The Super- 
Kamiokande detector response is given by R{Ee,E), 
which describes the probability that a recoil electron of 
energy E^ is reconstructed with energy E. The rate 
T = '"nooBc expected without oscillation is obtained by set- 
ting pe to 1. Since Pe is different for each neutrino species 
(pp, ^B, hep, etc.), the calculation must be repeated for 
each relevant species. Due to its recoil electron energy 
threshold of 5.0 MeV, Super-Kamiokande-I is only sen- 
sitive to ®B and hep neutrinos. The oscillation analysis 
was performed by two methods: the first method sub- 
divides the data sample in both recoil electron energy 
and solar zenith angle [Zenith Spectrum], while the sec- 
ond method uses only energy bins and searches for time 
variations by means of an unbinned likelihood [Unbinned 
Time Variation]. 



B. Analysis of the zenith spectrum 

If 3°"° {H""^) denote the expected *B (hep) neutrino 
induced rate in energy bin i and zenith angle bin z and 
Di_z the measured rate, the quantities 6, h and d are 
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FIG. 52: Dependence of Data/MC on solar zenith angle. The 
MC assumes a *B flux of 5.076 x 10^ /cm^s and a hep ly flux of 
32 X 10*' / cm^s. The six panels correspond to the energy ranges 
5.5-6.5 MeV, 6.5-8.0 MeV, 8.0-9.5 MeV, 9.5-11.5 MeV, 11.5- 
13.5 MeV, and 13.5-16.0 MeV. Overlaid is a typical LMA so- 
lution (red) and a LOW solution (green), where the assumed 
neutrino fluxes were flt to best describe the data. 



defined as follows: 



and 



D, 



Bi.z + Hi, 



For d, the rate calculations in the denominator use a 
full Monte Carlo simulation of the Super-Kamiokande 
detector. Figure [3] shows the average d in 21 and 8 
energy bins. In Figure the individual di^z are shown 
for 6 solar zenith angle bins. Now all zenith angle bins 
are combined into vectors. The rate difference vector 

aI(/3, t]) = (^P -1)1 + ri-~hi^ X f{E,, 63,63, 6r) - dl 

allows for arbitrary total neutrino fluxes through the 
free parameters P and 77. The combined rate predictions 
(3bi + rjhi are modified by the energy-shape factors 

fiE„SB,Ss,dB.) - fB[E„5B) X fs{E^,5s) X fBiE^Sn) 

with 6b describing the neutrino spectrum shape un- 
certainty, 5s describing the uncertainty of the SK energy 
scale (0.64%) and Sr describing the uncertainty of the 
SK energy resolution (2.5%). The n x n matrices Vi (n is 
the number of zenith angle bins) describe statistical and 
energy-uncorrelated systematic uncertainties. Those sys- 
tematic uncertainties are assumed to be fully correlated 
in zenith angle. 

To construct the Vi, we transform the rate vectors from 
the zenith angle basis to a new basis: the first compo- 
nent gives the average rate, the second component the 



day/night difference, the third component the difference 
between the rate of the mantle 1 bin and the other night 
bins, etc. This is achieved by multiplying by the non- 
singular matrix S which depends on the livetimes of each 
zenith angle bin: A/ = S" • 



-1 




^2/^1- 
L2I L2- 
-1 



V 







-1 



1 



J 



where is the livetime of zenith angle bin z and L^^n = 
J2^'=z ^z'- The quadratic form A^-^ • Vi~^ ■ Ai trans- 
forms as = S-V-S^ {{S-VrS^)-^ ^ S^-^-V-^-S-^). 
Since in this new basis, the first component of A/ is 
the livetime-averaged rate difference of energy bin i, 
the energy-uncorrelated systematic uncertainties af^ (as- 
sumed to be fully correlated in zenith angle) can simply 
be added to the statistical uncertainties S ■ gtat ' '■ 



S-V,-S^ = S- y,,stat • S'^ 



••• 



\ ••• 0/ 



Of course, Vi gtat is diagonal. If ai^u is asymmetric, two 
error matrices Vi+ and Vi- are constructed. The sign of 
A- J (livetime-averaged rate difference) decides, if A^^ • 
V~^} ■ A^ or A^^ • l/-!^ • A^ contributes to the 



i=l 



Vr 



(9.2) 



For any given parameters 5k, Xo is just a quadratic 
form in the neutrino flux factors [3 and 77 and can be 
written as 



xl 



Xn 



$ Co $ with 



(h — ( ^ /^min I 
\V-VmmJ 



where the 2x2 curvature matrix Co is the inverse of 
the covariance matrix for /? and rj. This matrix Co must 
therefore be inversely proportional to the combined un- 
certainty of all data bins: 



1 



with 



stat, 



stat, 



z = l 



1 



stat,; 



However, the energy-uncorrelated systematic uncertain- 
ties (Ti^u do not reflect the total uncertainty of the SK 
rate. For example, the uncertainty in the fiducial volume 
due to a systematic vertex shift cancels for the spectrum 
shape and is therefore not included in ai^u- If the spec- 
trum data is used to constrain the total rate, the total 
uncertainty of that rate (ctq) therefore neglects that part 
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FIG. 53: Predicted solar zenith angle variation assuming a 
typical LMA solution, Am^ = 6.3 x 10"'^ eV^ and tan^ 8 = 
0.52 (which is the best fit to the SK energy spectrum, rate, 
and day /night variation). 



(called (Tr) of the systematic uncertainty which cancels 
for the spectrum. While ar is fully correlated in both 
energy and zenith angle, it is not covered by the energy- 
correlated uncertainties either, which reflect only the un- 
certainties in the neutrino spectrum, the SK energy 
scale, and the SK energy resolution. To take ar into 
account Xo is modified to 



$ ^ Ci $ with Ci 



X Co 



which has the same minimum as Xo but allows an en- 
larged range of parameters rj. The total for the SK 
zenith spectrum shape is then 



XsK = Min( xliP,V,SB,6s,5R)+ 



(It) + (If ) + (If) ) 



(9,3) 



where all 6k as well as /3, 77 are minimized. To constrain 
the flux, the term is added. 



bin i, and r^^ the predicted time-averaged rate. Fig- 
ure shows the expected solar zenith angle dependence 
Ti (cos 6*2) in each energy bin for a typical LMA solu- 
tion. The recoil electron spectrum enters this likelihood 
through the weight factors Yi of Eq. IjS.ll) . so the likeli- 
hood needs to add terms like those in Eq. (|9.3|l to account 
for energy-correlated systematic uncertainties. The def- 
inition of energy bins is the same as in Eq. (|8.1|l : This 
analysis uses 21 energy bins compared to 8 energy bins 
for the zenith spectrum. 

Due to the large number of solar neutrino candidates, 
the maximization of this likelihood is too slow to be prac- 
tical. We therefore split it into 

log/: = (log/: - iog/:„.) -f iog/:„ 

where C^^ is the likelihood without the time variations, 
which means that it depends only on the average rates 
in each recoil electron energy bin. We can cast this in 
terms of a = — 21og/:: 

= Axi + xl with Axi = -2(log/: - log/:„) (9.5) 

The first term of Ax^.^ uses the time-dependent solar sig- 
nal factors of Eq. (|9.4|l while the second term is the same 
as Eq. H8.1|l . The spectrum weights Yi occurring in that 
equation are determined from the oscillated predicted 
spectrum. This predicted spectrum is formed using the 
^B and hep neutrino fluxes as well as the systematic un- 
certainty parameters 5b, a b and 5s- The values of these 
five parameters result from a fit to the (time-averaged) 
SK spectrum data. The dark-gray areas in Figure [5^ are 
excluded at 95% C.L. by this Ax?,. 

All systematic uncertainties are assumed to be fully 
correlated in zenith angle or any other time variation 
variable. Since the sensitivity is dominated by statistical 
uncertainties, this assumption is not a serious limitation. 
We found the statistical likelihood in each energy bin i 
to be very close to a simple Gaussian form, so xL takes 
the same form as Eq. (jnSIl but replacing Eq. ||n!21l with 



1=1 



(9.6) 



where A^ has only one zenith angle bin (between — 1 < 
cos 6*2 < 1) 



C. Unbinned time- variation analysis 

In a different approach to combine time-variation and 
spectrum constraints we form the total likelihood C mod- 
ifying the solar signal factors Sy of Eq. H8.ll) to 



A,(/5,77) = {P-h + T^-hi) X f{E„5B,5s,5R)^d, 

and of = erf ^^^^ + '''?« is the total energy- uncorrelated 
uncertainty in that bin. The regions colored in light-gray 
of Figure 15^ are excluded at 95% C.L. using xL without 
the term constraining the ®B neutrino flux. 



Sij = p{cos9ij,Ej) X 



(9.4) 



D. Oscillation constraints from SK 



where tj is the event time, ri{t) the predicted time de- 
pendence of the solar neutrino interaction rate in energy 



The combined excluded areas (from x'^ = + ^X?v) 
at 95% C.L. are shown inside the solid line in the left 
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FIG. 54: Left: excluded regions from the SK spectrum (red), 
SK rate time variations (green) and both combined (sohd hne) 
at 95% C.L.. Right: excluded regions from the SK zenith 
spectrum (green) compared to the SK spectrum combined 
with the time variation likelihood (red) at 95%C.L.. 



The best fit is in the LMA region at tan^ Q = 0.52 and 
Am2 = 6.3 X 10-5 eV^ (Am^ = 7.6 x lO'^ eV^) for the 
unbinned time- variation analysis (zenith spectrum anal- 
ysis) where a day/night asymmetry of -2. !%(-!. 5%) is 
expected and -1.8± 1.6% (-1.7± 1.6%) is found. (Here, 
the uncertainties do not include systematic effects) . The 

is 17.3 for 20 degrees of freedom (63% C.L.). The 
flux is fit to 4.91 X 10^/cm^s. The xL of the SK spectrum 
and rate is 18.5 for 20 degrees of freedom (55% C.L.). The 

analysis of the zenith spectrum gives a minimum of 
39.0 with 43 degrees of freedom (65% C.L.). The flux 
is fit to 4.86 X 10^/cm^s. Figure [S3 also shows the x^ as 
a function of tan^ B alone where a Am^ is chosen for each 
tan^ B to minimize x^ '■ SK data excludes small mixing at 
more than 3tT. SK data also disfavors Am^ > IQ-'^eV^ 
and 2 x lO'^eV^ < Am^ < 3 x IQ-^eV^ (see the plot 
in the right panel, where x^ is minimized with respect 
to tan^ B). A gain, the unbinned time-variation analysis 
has more stringent oscillation constraints and favors the 
LMA region more strongly. 
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FIG. 55: Allowed areas from the rate-constrained SK un- 
binned time-variation analysis (left) and the SK zenith spec- 
trum (right) at 95% C.L.. The graphs at the top (and right) 
show the difference as a function of tan^ 6 (Am^) alone 
where the Am^ (tan^ 6) is chosen to minimize x^- The lines 
inside the contours show which Am^ (tan^ 6') is chosen. 



E. Combined oscillation constraints from several 
experiments 

Stronger constraints on Am^ result from the combi- 
nation of the SK measurements with other solar neu- 
trino data. The 95% C.L. allowed region of Figure ICT a) 
include in addition to the SK data the SNO mea- 
surements of the charged-current interaction rate (and 
day/night asy mmetry) of solar electron neutrinos with 
deuterons zu\; only LMA solutions survive. Overlaid are 
the contours allowed by the KamLAND reactor neutrino 
spectrum [31 ■ Figure l^ b) shows the allowed contours 
of a combined fit to SK, SNO and KamLAND data. SK 
and SNO remove the ambiguities in Am^ and tan^ B of 
KamLAND and tighten the constraint on the mixing an- 
gle. When the charged-current rates measured by Home- 
stake, GALLEX, and SAGE are included as well, the 
allowed LMA solutions are further reduced. However, in 
this case the fit relies on the SSM predictions of the pp, 
pep, CNO and ^Be neutrino fluxes. 



F. Solar day/night effect 



panel of Figure |^ The right panel of the same fig- 
ure compares these areas to the excluded areas from the 
zenith spectrum analysis. Both contours are quite sim- 
ilar, however the unbinned time-variation analysis has 
more stringent limits on the time-variation (note the re- 
gion 10"^ eV^ < Ato2 < IQ-^ eV^ at tan^ 6 « lO^^ 
and tax? B > 1). Due to the larger number of energy 
bins, the likelihood analysis has also larger excluded ar- 
eas in the vacuum region and at Am^ « 10~^ eV^. If we 
include the last term in Eq. (|9.6|l . we get the allowed re- 
gions (95%C.L.) shown in Figure ICT which depend on 
the total ®B neutrino flux measurement of SNO [l^ . 



Even though KamLAND provides by far the best con- 
straint on the solar Am^ it is interesting to study the 
solar constraints as well. The upper bound arises from 
the ratios of the electron-neutrino elastic scattering rate 
in SK, the charged-current neutrino interactions in SNO 
and Homestake, and the neutral-current neutrino inter- 
actions in SNO. The lower bound is due to the solar neu- 
trino day/night effect. As described in we flt the 
amplitude of this variation to our data and compare it 
to the expected amplitude. For this fit, only the am- 
plitude of the day/night variation is varied, the shape is 
fixed to the calculated shape for a particular Am^ and 
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FIG. 56: (a) Allowed area at 95% C.L from the combination of SK and SNO. The graphs at the top (and right) show the 
difference as a function of tan^ 6 (Am^) only; the solid line is the SK/SNO fit, the dotted line results from KamLAND 'l&l 
data. The lines inside the contours show which Am^ (tan^ 6) is chosen, (b) Allowed area at 95% C.L. from the combination of 
SK/SNO (gray) and all solar data (dark gray), and from the combination of solar and KamLAND data (light gray). Note the 
linear scale of both axes. Only LMA-I remains allowed. Overlaid is the region allowed by KamLAND alone . 



tan^ 9 (see Figure ESI for a typical LMA solution) . De- 
manding consistency between expected and observed am- 
plitude, we obtain a range of 5 x 10~^ eV^ < Am^ < 
12 X 10~^ eV^ at the la level and a lower bound of 
3 X lO-'^ eV^ < Am^ at the 2a level (tan^ 6* = 0.44). 
This agrees very well with the KamLAND measurement 
of 7.9 X 10"^ eV^. Conversely, if we confine Am^ to 
the 3a range 7 x 10"^ eV^ < Am^ < 9 x 10"^ eV^ al- 
lowed by KamLAND, the amplitude fit varies only very 
slightly. In that range the measured SK day/night vari- 
ation amplitude corresponds to a day/night asymme- 
try of -1.7%±1.6%(stat)j:}:^(syst)±0.04%(Am2) while 
the expected asymmetry ranges from -1.7% to -1.0%. 
At the SK best-fit Am^ = 6.3 x 10"^ eV^, the 
day/night amplitude fit corresponds to the asymmetry 
-1.8%±1.6%(stat)ii-|(syst). 



X. CONCLUSION 

Solar neutrino measurements using neutrino-electron 
scattering in the Super-Kamiokande detector are de- 
scribed. We obtained 1496 effective days of data in the 
time period of May 31st 1996 through July 15th 2001. 
The analysis threshold was 6.5 MeV for the first 280 
days, and 5.0 MeV for the remaining 1216 days. The ob- 
served interaction rate corresponds to a *B solar electron- 
neutrino flux of 2.35 ± 0.02 ± 0.08 x lO^cm-^sec-i. We 
searched for periodic time variations of this rate and 
found only the expected seasonal variation caused by the 
eccentricity of the earth's orbit. The energy spectrum of 
the recoiling electron is consistent with an undistorted 
solar ®B neutrino spectrum. Based on these results the 
solar neutrino oscillation analysis imposes strong con- 
straints on the oscillation parameters, selecting large 
mixing and favoring the Large Mixing Angle solution. 
The combination with data from other experiments con- 
firms the Large Mixing Angle solution and further re- 
duces the uncertainty in the oscillation parameters. 
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FIG. 57; Day/Night amplitude fit depending on Am^. The 
fit depends only very weakly on the mixing angle which is 
assumed to be tan'^ 6 — 0.44. Overlaid (red line) is the ex- 
pected amplitude and (blue lines) the 2a boundaries of the 
amplitude fit. 
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